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INTRODUCTION 


General 


The Rocoz data reduction programs transform data from 
the Rocoz photometer to ozone number density and overburden 
as a function of altitude. Required auxiliary data are the 
altitude profile versus time and, for appropriate 
corrections to the ozone cross sections and scattering 
effects, air pressure and temperature profiles. Air 
temperature and density profiles may also be used to 
transform the ozone density versus geometric altitude to 
other units, such as to ozone partial pressure or mixing 
ratio versus pressure altitude. 

A set of seven programs are used to accomplish this. 
The programs and the functions of each are as follows: 

RADAR: Reads the 7-track BCD tape, converts from 
English to metric units, converts from seconds after 
launch to seconds after the hour, and writes real 
numbers on the 9-track tape. 

LISTRAD: Displays data from the "PASSl" radar tape 
generated by the Eclipse Computer at WFF. 

RAW OZONE: Displays selected records from the 
digitized telemetry tape to enable a preliminary 
assessment of the instrument performance and data 
qual ity. 

EDIT OZONE: For each rotation of the filter wheel, 
this program extracts a value for each filter, 
corresponding values of the compensation channel and 
of the uncompensated signal, converts telemetry data 
on battery voltage and instrument temperature to 
engineering units, and provides some statistics on 
instrument performance. 

MERGE: Obtains altitude information correspond ing to 
each filter observation; computes corresponding solar 
zenith angles 

SMOOTH: Derives smoothed photometer output at integer 
altitudes for each filter, with statistical measures 
of deviations. 

PROFILE: Converts the smoothed photometer data 
(intensity or intensity ratios) to ozone number 


1-1 



and air densities for evaluation of Rayleigh 
scattering effects. The program then converts the 
ozone densities and geometric altitude to other useful 
units: ozone partial pressure , mixing ratio, pressure 
altitudes, and, at intervals corresponding to standard 
pressure levels, SBUV pressure levels, and Umkehr 
layers . 


History 


The first of the series of Rocoz photometers was flown 
by A. Krueger in the 1950's; data reduction and analysis was 
with the assistance of a desk calculator. In 1977-8, D. U. 

Wright of GSFC and E. H. Shaffer of Systems and Applied 
Sciences Corporation (SASC) implemented the data reduction 
on the IBM 360/91 computers in the Science and Applications 
Computing Center (SACC) at GSFC. These programs were 
documented in the “Rocoz Automatic Data Processing User's 
Guide” prepared by SASC under Contract NAS5-24278. 

These programs were revised by K. Tewari of P & P 
Industries to accommodate changes in tape format, modified 
to reduce uncertainties, and adapted to special projects, 
part icularly , the use of the Rocoz instrument on various 
balloon platforms. 

This latest revision, by E. Reed and G. Batluck of GSFC 
and S. Cooke of Republic Management Systems is tne result of 
a line-by-line examination of the programs. Major revisions 
i nc 1 ude : 

The frames in each wheel rotation cycle are identified 
and handled in terms of the time tag rather than an 
index number. 

The noise from sources other than anal og-to-d i g ital 
conversion is evaluated. 

Some auxiliary parameters (e.g., earth radius) are 
computed more precisely. 

Numerous changes were made to improve efficiency and 
to add to clarity. 

Smoothing is accomplished by a least squares fit to t 
In S = A ♦ B Ir. (hs-hb). 

The Profile program has been comoletely rewritten. 
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The effective ozone absorption coefficient calculation 
can be readily modified to include the temperature 
dependence of the ozone absorption. 

A number of tests were added to the Profile program to 
evaluate the self-consistency of the data set. 

"Pedigree" records were added to the EDIT, RADAR, 
MERGE, SMOOTH, and PROFILE programs to enable 
inclusion of relevant auxiliary information in the 
archived data set. The Pedigree records of each 
previous step are copied to the output tape such that 
the Profile tape has a complete set of Pedigree 
records giving the entire history of the processing of 
the data. 

This document describes the result of this latest revision 
and includes details of the tape and data formats. 


Implementation 


These programs are written in Fortran-IV (Fortran-66) 
as implemented on the IBM-3081 computer in the Science and 
Application Computing Center (SACC) at GSFC. Tapes are 
identified through the FTIO set of subroutines rather than 
through the Job Control Language. Other system and 
SACC-suppl ied subroutines are: 

CMOVE for efficiently moving blocks of data within the 
programs 

BCD5 for converting BCD characters to EBCDIC 

SHFTL and SHFTR for extracting the integers which give 
time information. 


(Continued) 
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Data Requirements 


The following information is needed for the complete 
analysis of the data from a Rocoz flight, using the 
GSFC set of programs: 

1. Rocoz Data Sheet, identifying the vehicle and 
payload, giving statistics on time ">f flight and 
flight performance, and Dobson observations. 

2 Calibration data for the Rocoz photometer giving: 
Filter set number, when calibrated 
Calibration of filter, AO, AX, A2 , and/or filter 
shape 

- Zero offsets for SO, SI, S2, and S3 
Battery voltage calibration (word 6) 

Temperature sensor calibration (marker pulse 
height) 

3. Rocoz telemetry tape, 9-track 

4. Radar tape, 9-track, altitude, latitude, longitude, 
and time 

5. Data sonde information: 

Flight identification and launch time 

- Temperature vs altitude 

6. ECC data 

- Flight identification and time 

- Altitude, ozone partial pressure, ozone column, 
air pressure and air temperature versus time 

- Radar track of ECC sonde 

7. Strip chart of telemetry data (desireable, but not 

absolutely necessary) 
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RAW - Outline of Chapter 

Overview II-2 

Input Data 1 1 - 3 
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RAW PROGRAM 


OVERVIEW 


The sole purpose of the RAW Program Is to read and dump sample 
sets of the digital data from computer facilities at Wallops 
Space Flight Center to validate contents of the tape before 
proceeding with the subsequent processing steps. 

The sample sets may be chosen by record (Version 1 of Software- 
Input Card Data: Card 1, Variable 3) or by Seconds of Time 

(Version 2). 

The actual processing of these data takes place when the 'EDIT' 
Program is run. Details concerning the digital input data and 
the approach to the analysis of same is provided in the 'EDIT' 
Program section of this document. A description of the RAW 
Telemetry Tape Format and related pertinent information may be 
found in the section entitled 'TAPE FORMATS'. 


i 
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RAW PROGRAM 


INPUT DATA 


Columns 

Format 

Content 

Vari abl e 

Card 1 

1-6 

6A1 

Input Tape # 

TAPEIN 

9-10 

12 

Input File # 

NF 

13 

12 

Version # 

IVER 

16-19 

14 

Flight # 

NFLIT 

Card 2 

1-3 

13 

# of Samples 

ISAMP 

6-9 

14 

Beginning Rec # 

I BEG I N 

12-15 

14 

End Record # 

ISTOP 

17-20 

14 

Increment 

INCR 


Example: 

010 0001 0010 0050 

Print 10 sample sections of the first 10 records every 
50 records. 

Card 3 

1-6 F6.1 Start Seconds TIMES 

9-14 F6.1 End Seconds TIML2 

The first 2 cards are necessary and the 3rd card is optional, 
used only if output is chosen by seconds of time, rather than 
by Record # as in Version #1. Version #2 was never used 
after FI i ght # 290. 
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RADAR PROGRAM 


RADAR - Outline of Chapter 


Overview II-6 

Input Data II-7 
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RADAR PROGRAM 


OVERVIEW 

The Radar segment of software is intended to process the 
Radar tape as received from Wallops Space Flight Center only 
to the extent of reformatting the data such that it will be 
acceptable to the 'MERGE' Program which Is run after 'EDIT'. 

During the 'REFORMATTING', Radar data items not required in 
later processing were bypassed and a parity check was 
effected to properly locate the negative sign Intended for 
selected data words. The resulting file is written 
unformatted onto a tape using the FTIO (Fortran I/O) package. 

This 'REFORMATTING' effort was not necessary after Flight 
316 when we began receiving Radar tapes which could be used 
in the 'MERGE' processing phase in their original form. 
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RADAR PROGRAM 


INPUT DATA 


Columns 

Format 

Content 

Variable 

Card 1 

2-7 

6A1 

Input Tape 

TAPE I N 

14-15 

12 

Input Unit 

I'JNIT 

19-20 

12 

File # 

NF 

24-25 

12 

# Files 

NFILE 

30-36 

F7.2 

# Seconds after Hr 

STIME 

40-42 

13 

FI Ight 

I FLT 

46-47 

12 

Date-Year 

I DATE ( 1 ) 

49-50 

12 

Month 

IDATE (2 ) 

52-53 

12 

Day 

IDATE(3) 

Card 2 

2-7 

6A1 

Output Tape 

TAPOUT 

14-15 

12 

Output Unit 

IUNIT2 

19-20 

12 

Output File # 

NF2 

Card 3 

2-5 

F4.0 

Vehicle 

R(3) 

8-9 

12 

Launch Month 

LMO 

11-12 

12 

Launch Day 

EDY 

14-15 

12 

Launch Year 

LYR 

18-25 

F8. 1 

Launch Time (HHMMSS.S) 

R(5) 

28-34 

F7.3 

Latitude (Degrees 

R 6) 

37-44 

F8.3 

Longitude (Degrees) 

R(7) 

47-62 

4A4 

Launch Site 

SITE 


Vehicle Is either: 1. Rocket 

2. Balloon 


Launch Site is either: 01 Wallops Is., Va. 

02 Poker Flats, Alaska 

03 Natal , Brazil 

04 Marambio, Antarctica 

05 Palestine, Texas 


The above 

card input 

was 

used for flights up 

thru 

# 315. 

Starting 

with # 316, 

the 

following card input 

is used: 

Columns 

Format 


Content 


Variable 

2-4 

13 


Flight # 


IFLT 

7-13 

F7.0 


Launch Date (MMDDYY ) 


DATE 

16-21 

6A1 


Input Tape # 


TAPEIN 

24-25 

12 


Input File # 


NFILE 

Beginning 

with Flight 

# 

316, no radar output 

tapes 

are 

produced. 

Instead, the 

new version of radar 

tape 

is input 

di rectly 

into the Merge 

Program. 
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EDIT PROGRAM 


Edit - Outline of Chapter 


Overview 1 1 - 1 0 

Telemetry Data 1 1 - 1 0 

Identification of Cycles 1 1 - 1 1 

Selection of Filter Values 11-12 

Averaged Values 11-13 

Common Areas 11-14 

Input Data 11-15 
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EDIT PROGRAM 


Overview 


The primary function of the EDIT Program is to extract for 
each rotation (cycle) of the filter wheel, a value for each 
filter and the corresponding compensation word. In 
addition, the battery voltage monitor and instrument 
temperature sensor are read and averaged, the length of each 
rotation (cycle) is noted, and some statistics on noise are 
generated . 

The filter wheel contains four filters, designated SO, SI, 
S2, and S3, with the S3 filter at the shortest wavelength, 
SO at the longest. When each filter is fully in the optical 
path of the instrument, the output of the photodiode is 
sampled and held, and a constant value appears on the 
telemetry word for the remainder of the cycle. A marker 
pulse on a separate telemetry word (channel) indicates the 
beginning of each rotation, defined as the time 
corresponding to an observation through the SO filter. 


Telemetry Data 


Data is transmitted at rate of 8K bits/second, or 88.9 
frames per second where each frame consists of a 
synchron i zat i on word followed by 8 data words, each 10 bits 
long. A typical rotation requires about 650 ms, during 
which the data are sampled about 58 times. The telemetry 
words are defined as follows: 

Word Function 

1 Signal from the SO filter 

2 Signal from the SI filter 

3 Signal from the compensation detector 

4 Signal from the S2 filter 

5 Signal from the S3 filter 

6 Battery voltage monitor 

7 Marker pulse, temperature monitor 

8 Signal from detector, before compensation 

During the initial data reduction process each frame is 
tagged with time information (day, hour, minutes, seconds, 
and milliseconds) and grouped into records, each of which 
contains 110 frames of data. 

In this program, the input data is read in Subroutine BLOCK. 
The Header r>cord and Data Description records describing 
the tape are converted and displayed. The data words are 
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EDIT PROGRAM 


converted to IBM format by switching bytes. The three time 
words appended to each frame of data are broken apart and 
interpreted in terms of the date and time. 


Identification of Cycles 


The identification of the start of each cycle for each of 
the four filter words is done in Subroutine CYCLE. CYCLE 
transfers data from the input telemetry tape records to the 
array, KEEP, always maintaining 200 frames of data in KEEP, 
and calling Subroutine BLOCK whenever another telemetry 
record is required. Before moving a frame into KEEP, the 
data is tested to be sure that the times are increasing and 
that word 7 (marker pulse) channel is either near the 
expected base level or the expected marker level, thus 
identifying and discarding frames that show evidence of 
being noisy. The initial expected values (MPEAK and LPEAK) 
are provided on cards as obtained from inspection of a 
listing of the telemetry data, and are updated by a running 
average as the program processes the tape. The timing of 
the changes is diagrammed as follows: 

M 


SO- 

SO 

SI- 

51 

S 2 - 

52 

53 - 
S3 



The beginning of a set of values for word 1 (SO) is taken as 
the first frame after a base value and containing a marker 
value. The beginning of the word 2 set (SI) is the opposite 
transition, that is the first frame with a base value after 
a frame with a marker value. The end of the SO set of 
values is the frame with the last base value before a marker 
value. The time interval between the start of the SI set 
and the end of the SO set is divided into thirds; the S2 
set of words is assumed to start at the end of the first 
third; the S3 set of words is assumed to start at the end 
of the second third. The length of the SI set is taken to 
be 94% of the length of the SO set, S2 set is 88%, and S3 
set is 82% of the SO set. (The shorter times for some 
filters allow for uncertainties in the start times due to 
noisy data and for variations in the cycle time.) Thus, 40 
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to 60 samples are available for each of the four filters, 
unless diminished by noise and telemetry drop-outs. 


Selection of Filter Values 


Since the SO - S3 words are presented to telemetry by 
sample-and-hold circuitry, each channel should remain 
constant for a complete rotation. Inspection reveals that 
this is not the case: the first word is sometimes at the 

old, rather than the new level; the digitization error is 
typically ♦ one bit; noise either in the instrument and the 
telemetry system can introduce greater variability; 
drop-outs in the received telemetry signal can reduce the 
number of available samples. Subroutine AVE selects a value 
representative for the array of words available for each 
filter. 

The rules used by Subroutine AVE are: 

- If only one sample is available, use it. 

- If two samples are available and agree within 6 counts, 
take their average. 

- If three or more samples are available, a subset is chosen 
in which all the words are +3 counts of a test value. If 
this subset represents more than half of the samples 
available, the average of the values in this subset is 
used. If the subset represents less than half of the 
values, a new test value is chosen (the next word) and a 
new subset tested. 

- If the above rules do not result in an acceptable value, 
the word is set to -99., indicating that there is not an 
acceptable value for this filter during this rotation. 

Subroutine AVE is used not only to obtain representative 
values for each filter but also for the marker pulse height 
and for the marker pulse base level. 

For each of the four filters, Subroutine COMPWD finds the 
value of the compensation channel (word 3) corresponding to 
the time of the first sample for the filter. It chooses the 
frame near the start frame chosen in Subroutine CYCLE, for 
which word 8 (the uncompensated signal through the filter 
wheel) is at a maximum; word 3 of this frame is then chosen 
as representing the value of the compensation word at the 
moment that the filter wheel signal was sampled and held. 
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Averaged Values 


Running averages of the most recent ten values are 
maintained for the height of the marker pulse, the base 
level of the marker, and the length of the cycle and used 
for the detection and deletion of aberrant data. To do 
this. Subroutine Cycle calls Subroutine TENAVE as it nears 
the completion of processing for each cycle. 

The average battery voltage is obtained from a simple 
average of all the values of word 6 in the most recent 
record of telemetry data (110 frames) transferred into the 
array KEEP by Subroutine CYCLE. The instrument temperature 
is obtained from the value of the marker pulse in the cycle 
under consideration. The average value of all the 
compensation words in the cycle is obtained in Subroutine 
CMPWD . 

The average values for the temperature and the battery 
voltage are converted into engineering units in Subroutine 
CONVRT. Calibration data are obtained from the preflight 
testing and given to the program through DATA5 cards. The 
temperature data are read in at three temperatures and a 
standard temperature versus counts curve is adjusted to fit 
the calibrations at the three temperatures . This adjusted 
curve is then used to convert the counts to temperatures . 
Calibration values for the battery word are read in at two 
points, and a straight line relationship is assumed. 
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Common Areas 


Four labeled common areas are used: LABELI transfers 
initialization data from MAIN to Subroutine CYCLE; LABEL3 
contains the telemetry data, transferring it from BLOCK to 
CYCLE; LABEL5 contains the output data generated in various 
subroutines;, and CAL contains calibration data, making it 
available in MAIN for use in the Pedigree records and in 
Subroutine CONVRT for conversion of telemetry data to 
engineering units. 


COMMON/ LABEL 1 / IREC , DE LTAT ,MPEAK ,LPEAK , LI MIT 


Variable 

urn: — 


OELTAT 

MPEAK 

LPEAK 

LIMIT 


w 


R*4 

1*4 

1*4 

1*4 


Content 

Number of the last record to be 
processed (Statement is in Block, 
usually commented out.) 

Initial value for length of cycle, ms 
Initial value for marker pulse peak. 
Initial value for marker base level. 
Delta for acceptable marker values. 


COMMON/ LABE L3 / DATA ( 1 210) .INPUT (1220) 

DATA(1210) R*4 Input record in IBM format, no header. 

INPUT ( 1220)R*4 Input record from telemetry tape. 


COMMON/ LABEL5/ OUTX ( 20 ) 

OUTX ( 20) R*4 Output array. 


C0MM0N/CAL/B1 ,V1,B1,V2,T(3) ,VT(3) 


B 1 

R*4 

Battery voltage, first 

calibration pt. 

VI 

R*4 

Counts corresponding to 

Bl. 

B2 

R*4 

Battery voltage, second 

cal ibrat ion . 

V2 

R*4 

Counts correspond i ng to 

B2. 

T ( 3 ) 

R*4 

Temperatures (deg. C) for calibration. 

VT( 3) 

R*4 

Counts corresponding to 

T(3). 
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Input Data 


Card 

Column 

Format 

Content 

Variable 

1 

1-6 

6A1 

Location of input tape 

TAPEIN 


9-14 

6A1 

"INPUT" 

10 1 N 


17-18 

12 

Input unit 

IUNIT1 


21-22 

12 

Input file 

NF1 

2 

1-6 

6A1 

Location of output tape 

TAP0UT 


9-14 

6A1 

"OUTPUT" 

I00UT 


17-18 

12 

Output unit 

IUN. "2 


21-22 

12 

Output file 

NF2 


24-26 

13 

Length of output record, 80 

LEN2 

3 

1-5 

F6.3 

Initial length of cycle, sec. 

DELTAT 


8-10 

13 

Number of first record to be 
processed . 

IREC 


13-15 

13 

Initial value of marker pulse, 
counts. 

MPEAK 


18-20 

13 

Initial value of marker base, 
counts . 

LPEAK 


23-25 

13 

Acceptable range of marker. 

LIMIT 

4 

1-5 

F5.1 

Battery voltage. 

B1 


6-13 

F8.2 

Corresponding counts. 

VI 


14-18 

F5.1 

Battery voltage. 

B2 


19-26 

F8.2 

Corresponding counts. 

V 2 

5 

1-5 

F5.1 

Low temperature, degrees C. 

T(l) 


6-13 

F8.2 

Corresponding counts. 

VT ( 1 ) 


14-18 

F5.1 

Mid temperature. 

T(2) 


19-26 

F8.2 

Corresponding counts. 

VT ( 2 ) 


27-31 

F5.1 

High temperature 

T( 3) 


32-29 

F8.2 

Corresponding counts. 

VT ( 3 ) 

6 

1-5 

F5.1 

Flight number. 

0UTX ( 2) 


8-11 

F4.0 

Payload number. 

0UTX ( 3) 


14-15 

12 

Month Telemetry Tape received. 

IDATE(l) 


17-18 

12 

Day 

IDATE (2) 


20-21 

12 

Year 

IDATE (3) 
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Merge - Outline of Chapter 
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MERGE PROGRAM 


Overview 


To compute the ozone concentration it is necessary to have 
the signal from each filter as a function of altitude. In 
the Merge Program, the edited telemetry data and the radar 
data are merged on the basis of time (in seconds after t!ie 
hour) as the common variable. The solar zenith angle is 
computed and included in each record of the Merge tape. For 
the lower part of the flight, a table of counts fo r each 
filter as a function of temperature is constructed and 
displayed. The basic structure of the SASC/P&P version of 
the program is retained, with the addition of the Pedigree 
records, provision for direct use of the "PASS1" radar tape 
produced by Code 961 on the Eclipse computer at WFF, 
determination of counts versus temperature, and other minor, 
mostly cosmetic, changes. 

After mounting and reading the first records on the Edited 
Ozone Data Tape and the radar tape, the agreement in flight 
numbers is verified, and up to 3000 records of radar data 
are read, storing the time and location data in arrays. The 
solar zenith angle is computed every 60 seconds and stored 
in another array. Finally, ozone data are read, the 
altitude for each filter word interpolated from the radar 
arrays, the solar zenith angle interpolated for the time of 
this record, and the Merge records written on the output 
tape . 


Radar Data 


For flights through 1982 (through Rocoz flight 315) radar 
w ata was provided in the form of a "MESUP" tape, with 
smoothed data for location and velocities given in English 
units and written in character format on a 7-track tape. 
This is then converted to metric data in binary files on a 
9-track tape by the Radar Program. The resulting 
"Reformatted Radar Tape" is then suitable as input to the 
Merge Program. 

For flights beginning with Rocoz no. 316 (August 1983), 
radar data Is supplied by CODE 961 at WFF, using the Eclipse 
computer. This "PASS1" tape gives the location as a 
function of time, in meters and in degrees of latitude and 
longitude. The Merge Program has been modified to permit it 
to read this tape directly. 
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Both types of tapes from WFF provide time in terms of 
seconds since launch. These must be converted to seconds 
since the hour prior to merging with the telemetry data. 
For MESUP tapes this is done in the Radar Program; for the 
PASS1 tapes, this is done in the Merge Program. The Merge 
Program also creates a Pedigree record for the PASS1 tape. 


Zenith Angle Computation 


In Subroutine ZENITH the solar zenith angle with respect to 
the Rocoz payload is computed. The apparent right ascension 
and declination of the sun and sidereal time is provided 
through DATA5 cards; the location of the payload is 
obtained from the radar tape. If the radar data is from the 
MESUP tape, then the location in distance from the launch 
site must be converted to latitude and longitude. The earth 
radius required for this conversion is obtained from the 
latitude of the launch site by the formula (from Basic 
Physics of the Solar System, Blanco and McCuskey, 
Addison-Wesley Publ., 1961, p. 83): 

R0-6378388. ( 1-3 . 367x 10 ' 3 ( s i n ( 1 at ) ) 2 +7 . 1 xlO " 6 ( si n ( 2 1 at ) ) 2 ) 

where R0 is the radius of the earth at the launch site and 
lat is the latitude of the launch site. The location of the 
payload is then: 

latp * latO ♦ 57.29578 rns/RO, and 

lonp * lonO + 57.29578 rew/(R0 cos(latO) ) 

where latp and lonp are the latitude and longitude of the 
payload, latO and lonO are the latitude and longitude of the 
launch site, and rns and rew are the distance (meters) north 
and east of the launch site. With the appropriate 
trigonometry, the solar zenith angle as a function of 
elapsed time and payload location can then be computed. See 
the program listing for the steps in detail. 


The Merge Process 


After the radar is stored in arrays, and the solar zenith 
angle computed for each minute, the Edited ozone data is 
read In, a record at a time. Provision is included for 
skipping “bad" sections of data as defined by the input 
DATA5 cards, although in recent flights this has not been 
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used. The altitude at the time of each word is interpolated 
from the radar data and the solar zenith angle at the time 
of the SO filter is interpolated from the array of zenith 
angles and are placed in an output record for the Merge 
tape. The interpolation is done by Subroutine AITKEN, which 
uses a four-point scheme. Also on the Merge tape are the 
time for the SO filter, the value for each filter during 
this cycle, the corresponding compensation word values, and 
information concerning the operation of the photometer, 
namely, the length of the cycle, the number of telemetry 
samples (frames) used in evaluating the cycle, the 
temperature of the photometer, and the battery voltage. 
Processing is terminated when the end of either the radar 
data or the Edited ozone data is reached. All values 
written on the tape are also listed for inspection prior to 
running the Smooth Program. 


Counts vs Temperature 


To assist in the determination of the zero offset values for 
filters S2 and S3 as a function of temperature. Subroutine 
TMPAVE is called for each Merge record as it is generated. 
For altitudes less than 45 km, S3 less than 100 counts, and 
temperatures between 39.74 and -10 degrees C., the counts 
for S3, S2, and SI are accumulated in 0.5 degree buckets. 
At the conclusion of the Merge Program, Entry PRINTM of 
TMPAVE is called to calculate the averages and print a table 
of the results. In calculating the average count, the 
maximum and minimum values available in each bucket are 
discarded . 
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Common Areas 


Only one common area is used, named Z, to transfer the 
astronomical data, time, and payload location data from MAIN 
to Subroutine Zenith where they are used in the computation 
of the solar zenith angle. The resulting solar zenith angle 
is returned to MAIN via this same common area. 


Variable 

ll21 

Content 

RH ,RM ,RS 

R*4 

Apparent Right Ascension in hours, 
minutes, and seconds 

RACOR 

R*4 

Change of Apparent Right Ascension, 
seconds per day 

DD ,DM , DS 

R*4 

Apparent Declination in degrees, 
minutes, and seconds. 

DECOR 

R*4 

Change of Apparent Declination, 
seconds per day. 

SH.SM.SS 

R*4 

Apparent Sidereal Time in hours, 
minutes, and seconds. 

GMT 

R*4 

Greenwich mean time in hours. 

REW.RNS 

R*4 

Distances east and north from launch 
site in meters; if M0DE*2, degrees 

SOLARZ 

R*4 

Solar zenith angle, degrees 


Input Data 


The Merge Program uses either the Reformatted Radar Tape or 
the PASS1 Radar Tape and the Edited Ozone Data Tape as 
input; detailed descriptions of these tapes are given in 
the section on Tape Formats. 0ATA5 "cards" are used to 
Identify the input and output tapes and files, to state the 
flight number, and other data needed to construct a Radar 
Pedigree record, to provide the astronomical data, and to 
delineate "bad" time periods on the Edited Ozone Data Tape. 
The Apparent Right Ascension and Declination of the Sun may 
be obtained from tables for Oh Ephemeris Time In Section C 
of the Astronomical Almanac and the Sidereal time from 
tables of Universal and Sidereal Times in Section B. 
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Card 

Columns 

Format 

Content 

Variable 

1 

2-7 

6A1 

Edited 

Ozone Data Tape ID 

EDITED 


12-16 

15 

H 

" file number 

NF 


17-21 

15 

II 

“ unit number 

I UNI T 

2 

2-7 

6A1 

Merge 

Data Tape ID 

MERGED 


12-16 

15 

H 

" file number 

NF1 


17-21 

15 

H 

" unit number 

IUNIT1 

3 

2-7 

6A1 

Radar 

Data Tape ID 

RADAR 


12-16 

15 

ii 

" file number 

NF2 


17-21 

15 

ii 

" unit number 

IUNIT2 


25-26 

12 

Mode 

reformatted radar 

MODE 





«2 PASS1 radar 


4 

1-5 

F5.1 

Rocoz 

flight number 

PED ( 2 ) 

5 

1 

11 

Locat ’ 

on 

L0CA 




«1 

Wallops Island 





*2 

Church 1 1 1 





-3 

Natal 





*4 

Primrose Lake 





*5 

Poker Flats 





«6 

Marambi o 



28-33 

F6.1 

Time correction for Edited 

RESTR 




Ozone, usually 0., seconds 



36-38 

13 

Begin 

Edited Ozone record no. 

IBEG 


41-42 

12 

Number 

of bad sections of 03 data NBAD 




maximum of 20 



43-72 

30A 1 

Rocoz 

no . , date , and time 

LAUNCH 

6 

11-20 

F10.2 

Apparent Right Ascension, hour 

RH 


21-30 

F10.2 

H 

M " minutes 

RM 


31-40 

F10.2 

II 

" " seconds 

RS 


41-50 

F10 .2 

II 

Correction sec/day 

RAC0R 

7 

11-20 

F10.2 

Apparent Declination, degrees 

DD 


21-30 

F10.2 

H 

" minutes 

DM 


31-40 

F10.2 

M 

M seconds 

DS 


41-50 

F10.2 

N 

Correction sec/day 

DECOR 

8 

11-20 

F10.2 

Sidereal Time hours 

SH 


21-30 

F10.2 

M 

M minutes 

SM 


31-40 

F10.2 

N 

" seconds 

SS 


41-50 

F10.2 

Launch 

time, hour, GMT 

HOUR 




from which seconds are counted 


9-28 

1*6 

F6.1 

Start 

of section to be skipped 

BAD1 




seconds after hour 



9-14 

F6.1 

End of 

section to be skipped. 

sec BAD2 
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Overv i ew 


The Smooth Program obtains a value representing the flux 
through each of the filters at each integer km of altitude. 
These values are obtained from lines which are a least 
squares fit to the data in the vicinity of each altitude. 
If requested. Smooth will do a similar computation for the 
ratio of flux to that through the SO filter. The results, 
with related information, are written on the Smooth Tape for 
input to the Profile Program and also presented in the form 
of tables and plots for examination and evaluation. 

For input, the Smooth Program requires the Merge data tape, 
values for the zero offsets and their temperature 
dependence, the minimum acceptable value for the 
compensation word, and the desired altitude range. 

The Smooth Program descritjd here is based on the code 
supplied by P & P Industries; much of the descriptive 
material in the SASC "Rocoz Automatic Data Processing User's 
Guide" of January 1979 is still applicable. The principal 
changes are fitting to a logarithmic function instead of a 
polynomial and basin? the error on the fit of the data 
rather than assuming a fixed count. Additions include a 
test of the compensation word and the use of the temperature 
coefficient of the zero offset. 

MAIN reads the input from DATA5 and from tape, places all of 
the Merge data on a disk, calls SELECT to identify and 
process the data in the vicinity of each altitude level, 
calls LOGFIT to obtain values representative of each 
altitude level, calls OUTPUT to write the output in tabular 
form and on a Merge tape, and finally calls GRAPH to plot 
altitude and intensity. 

All tape and disk input and output are through calls to the 
FTI 0 package available on the IBM-3081 computer at GSFC 
Science and Applications Computing Center. Graphs are 
prepared with the assistance of the WOLFPLOT package 
likewise available on that computer. 
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Selection and Preparation of Data 


Initially, the data records are copied from the Merge tape 
to a disk data set. The data are assumed to be ordered as 
descending altitude. (For analysis of data from the ascent 
portion of a balloon flight, this program would require 
modification.) MAIN chooses integer altitudes, beginning at 
the top of the altitude range read from DATA5, and calls 
Subroutine SELECT. 

SELECT searches the disk data set for Merge records in the 
vicinity of this center altitude. Equal numbers of records 
above and below that for the center altitude are chosen 
except as limited by the availability of data. Collection 
of records continues until 100 records are available and the 
altitude interval is at least two km, or until a maximum of 
800 records are collected. The average values for cycle 
length and the compensation word are computed from the 
values given in the 99 Merge records centered at the 
selected altitude. If requested, records are created for 
which the fluxes for the S3, S2, and SI filters have been 
divided by the SO flux. 

Not all of the data are used. Excluded are: records more 
than 0.5 km above the top of the requested range; filter 
words with a value of less than 2.0 counts (before zero 
correction); and filter words for which the corresponding 
compensation word is below the minimum acceptable value. 

Subroutine SELECT corrects all the remaining filter words 
for zero offset. The zero offsets for the SO and SI filters 
are taken from the records of the laboratory preparation of 
the payloads, and are usually 0.0. Values for the S2 and S3 
filters are obtained from inspection of the Merge listing of 
the S2 and S3 values at low altitudes. A plot of these data 
versus temperature usually reveals a non-zero temperature 
coefficient. Values at two temperatures are placed in 
DATA5 ; In SELECT, a temperature dependent offset is computed 
for each filter word. 

Finally, In Subroutine L0GFIT, filter words which are more 
than two standard deviations from the initial fitted line 
are discarded, and the constants for the line recomputed. 
The number of records selected and the number of filter 
words used are included in the Smooth output record. 
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Smoothing and Error Estimates 


In Subroutine LOGFIT (based on a program given by 
P.R.Bevington in Data Reduction and Error Analysis for the 
Physical Sciences, Mc6raw-Hill, 1969, p. 104-105, 180-185), 
the arrays of flux (or ratio to SO) and altitude are fitted 
to a line of the form 

In I ■ In A ♦ B (h • hb) 
where I is the flux or flux ratio, 

A is the value of I corresponding to hb) 

h is the altitude, and 

hb is the lowest altitude in this array. 

Normally each value is given equal weight (MODE = 0) , 
although the code includes provision for other weighting 
schemes. Values for A, B, and the standard deviation for 
the array are calculated, and all points which are more than 
2 standard deviations from the calculated value of In I are 
discarded. The remaining values in the array are used in a 
recomputation of A, B, correspond i ng standard deviations 
(SIGMAA and SIGMAB), and the correlation coefficient. 

This treatment is not appropriate for all arrays of data, 
particularly those at large optical depths for which the 
signal is essentially constant with altitude. Such arrays 
are detected by testing DELTA and variance for near-zero 
values. They are evaluated by obtaining the average value 
of I over the altitude range and calculating the 
corresponding standard deviation. Here, also, values of I 
which are more than 2 standard deviations away from the 
average are discarded and the average recomputed. When this 
method is used, the values of A, B, SIGMAA, SIGMAB, and the 
correlation coefficient in the Smooth Data Record are set to 
zero . 

Subroutine SHAPE is provided which may be called from 
Subroutine L0G r IT (see the comments immediately after ISN 
0053) to list the residuals for all the points of an array. 

To provide a measure of the extent of smoothing, the number 
of points used and the upper and lower altitude limits of 
each array is included in the Smooth Data Record. 

In an earlier version of this program, provision was 
included for associating an error with each point and using 
this error in the evaluation of the error in the fitted 
line. This feature is no longer used, although the error 
arrays are still present. 
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Graphical Displays 


Subroutine GRAPH, called from MAIN, uses the WOLFPLOT 
package to prepare plots of the smoothed value of the signal 
(or ratio to the SO signal) as a function of altitude at 
each km, using values from the SAVE array. These are 
displayed as a printer plot and can also be written to tape 
for use with higher resolution plotters. 


Common Areas 


Two named common areas are used, one for the arrays of data 
associated with a particular altitude level, and the other 
for the output records. 

COMMON/ L ABE LI /N,S( 4, 800) ,H ( 4 ,800 ) ,ER ( 4 ,800 ) 

Found in MAIN and Subroutine SELECT. 


N 

1*4 

Number of points 
interval 

in this fitting 


S ( 4 , 800 ) 

R*4 

Number of counts (telemetry) for 
of the four filters. 

each 

H ( 4 ,800 ) 

R*4 

Altitude for each 

of the points 

(km) 

ER ( 4 , 800 ) 

R*4 

Error associated 

with each point 

. 


COMMON / LABE 12/ SAVE ( 20 ,50,4) 

Found in MAIN and Subroutine OUTPUT. 

SAVE ( 20 ,50 ,4 ) R*4 Array containing all output records. 


Input Data 


The Smooth Program uses the Merge Tape as the source of data 
from the Rocoz photometer as a function of altitude; the 
detailed description of this tape is given in the section on 
Tape Formats. DATA5 "cards" are used to identify input and 
output tapes and files, run options, altitude range, minimum 
compensation channel value, and zero offsets to be applied 
to each filter. The data in terms of telemetry counts is 
always evaluated; upon request, the data as the ratio to 
the SO signal is also smoothed. 
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The altitude range is obtained from inspection of t.ie 
listing from the Merge program. The upper limit is the 
altitude at which the attitude is relatively stable, usually 
several km lower than the ejection altitude; the lower 
limit is the lowest integer altitude for which data are 
available. The minimum compensation word value is usually 
set equal to one-half the minimum compensation value 
observed, such that data more than about 60 degrees from the 
Rocoz optical axis can be discarded. Zero offsets for SO 
and SI are obtained from laboratory measurements; for S2 
and S3, the data at the lower altitudes are plotted as a 
function of temperature, a straight line drawn, and values 
read at two temperatures. 


Card Columns Format Content 


Variable 


1 

1-6 

6A1 

Input Tape Identification 

TAPEIN 


9-10 

12 

Input Tape File number 

NF 


13-14 

12 

Input Tape Unit 

IUNIT 

2 

1-6 

6A1 

Output Tape Identification 

TAP0UT 


9-10 

12 

Output Tape File number 

NF2 


13-14 

12 

Output Tape Unit 

IUNIT2 

3 

1-2 

12 

Type of weighting, usually 00 

MODE 


45 

11 

=0: no ratios; -1 divide by SO 

LRAT 

4 

1-6 

F6.1 

Zero offset (counts) for SO 

OF F ( 1 ) 


9-14 

F6.1 

Zero offset (counts) for SI 

OFF ( 2) 


17-22 

F6.1 

Zero offset (counts) for S2 

OF F ( 3 ) 


25-30 

F6.1 

Zero offset (counts) for S3 

OF F ( 4 ) 

5 

1-6 

F6.1 

Temperature (deg. C) for above 

TA 


9-14 

F6.1 

Temperature for next two values 

TB 


17-22 

F6.1 

Zero offset for S2 at 2nd temp. 

S2B 


25-30 

F6 . 1 

Zero offset for S3 at 2nd temp. 

S3B 

6 

1-4 

F4.0 

Rocoz flight number 

PED ( 2 ) 
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Overview 


The Profile Program obtains data for the Rocoz photometer 
from the Smooth tape, reads selected values for the 
Oatasonde, the ECC sondes, and the Dobson spectrophotometer, 
and reads filter calibration values and related tables. 
The program uses these data in the computation of ozone 
profiles for each filter and also creates a composite 
profile based on data from all the filters. With the air 
temperature and pressure data available from the Oatasonde 
and the ECC sonde, the ozone profile is then be expressed in 
units such as partial pressure and mixing ratios. Various 
aspects of the data set are examined for self-consistency. 
Finally, the basic data and the derived profiles are written 
on a Profile tape, intended as the archival product of the 
flight. The Profile Program does not create plots of the 
data but does include a number of printed tables. These 
functions are distributed among the 21 Subroutines as 
f ol 1 ows : 

(1) Assembly of the Rocoz data and correlative 
observations with the required calibration data. 
(CARDS, TAPE, INARRY, SHAPE) 

(2) Computation of related data. 

(MODEL, DATSND , SECANG) 

(3) Computation of ozone profiles and composite 
profile. 

( OZDENS , SLANTA, ALFEFF, COMPST) 

(4) Evaluation of Rocoz data set. 

(CHECKS, OVERLP, LINFIT) 

(5) Computation of related parameters and display of 
observations in tables. 

( RELATD , PRTWRT) 

(6) Placement of the ozone profiles and related data 
on the Profile Tape. 

( TAPWRT) 

(7) Display of various arrays as desired for test and 
diagnostic purposes. 

( PRTPRF , PRTPED , CFTWRT , PRTSAV) 

The general structure is shown in the Block Diagram of the 
Profile Program. 
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CARDS 

TAPt 

MODEL 

DATSND 

SECANG 

OZDENS 

SLANTA 

ALFEFF 

INARRY 

SHAPE 

COMPST 

CHECKS 


Brief Description of Each Subroutine 


reads set-up information (e.g. altitude range for 
each filter) and correlative data from the 
Datasonde, the Dobson spectrophotometer , and the 
ECC sondes. 

reads Rocoz data from the Smooth tape. 

contains an ozone model, with derived overburdens. 
The present Krueger-Minzner mid-latitude model 
should be replaced with one dependent on latitude 
and season. 

calculates densities and pressures correspond i ng to 
the datasonde temperature profile. 

obtains the secant of the solar zenith angle for 
each altitude level; for angles over 60 degrees, 
the Chapman function is computed instead. 

computes ozone density and column content for each 
altitude level for each filter, based either on 
observed fluxes or their ratio to the SO flux. 

computes the slant air mass (atmospheres) and the 
slant ozone overburden (atm-cm) at the three 
altitude levels under consideration in OZDENS. 

computes effective ozone absorption and Rayleigh 
scattering coefficients, given filter transmission, 
temperature, slant overburdens, and ozone absorption 
coefficients. 

obtains solar flux, diffuser transmissions, detector 
response, Rayleigh scattering coefficients, and 
ozone absorption coefficients from external data 
sets, determines range of wavelengths for which data 
are available, and sets up these arrays for use in 
computing the effective absorption coefficient. 

obtains the measured filter transmission. 

obtains a weighted (inversely to associated error) 
of the available Rocoz ozone data at each altitude 
level. ECC data arc not included. 

tests the internal consistency of the data related 
to this profile by (1) comparing ozone values at a 
given level from different filters, (2) computing 
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OVERLP 

UNFIT 

RELATD 

PRTPRF 

CFTWRT 

PRTPED 

PRTSAV 

PRTWRT 

TAPWRT 


F 7 


ozone overburden from the Composite Profile as a 
function of altitude and comparing it with the 
Model, (3) comparing the measured ozone densities 
with those from the Model, (4) comparing the 
(Rocoz+ECC) column content with that measured by the 
Dobson instrument, (5) comparing various estimates 
of the rate of change of In I, (6) comparing the 
variation of the intensity in the compensation 
channel with that expected due to changes in the 
solar zenith angle and scattering, and (7) computing 
the drag coefficient as a function of altitude 
(large variations may indicate a defective Starute 
or errors in the radar altitude data). 

called by CHECKS to assemble and analyze the ozone 
densities at each altitude level. 

called by OVERLP to compute a straight line fit to 
the data at each level . 

incorporates air temperatures and pressures, and 
computes ozone values within standard layers and at 
standard pressure levels. 

displays the contents of the COMMON/ PROFLS/ which 
includes the Datasonde, the ECC, and the model 
prof i 1 es . 

displays the contents of the arrays generated by 
INARRY and FSHAPE , used ir. the computation of 
effective absorption coefficients. 

displays the contents of the Pedigree array. 

displays the contents of the SAVE array for each 
filter, listing the ozone profile and values used in 
its derivation. 

formats arrays for hard copy printout with 
appropriate labels and column headings, including 
the arrays generated by RELATD. 

writes the data selected for archiving on a tape. 
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Observations and Supplied Values 


In the analysis of observations from a Rocoz photometer, 
additional data are used: various constants, an air 
temperature profile from a Datasonde, air pressure and ozone 
profiles from an ECC sonde, total ozone from a Dobson 
spectrophotometer, and a model ozone profile. The flight 
plans require that a Datasonde and two ECC sondes be 
launched close to the time of the Rocoz launch, that all 
sondes be tracked by radar, and that the Dobson 
spectrophotometer be operated according to its usual 
schedule. The raw data from the Datasonde, the ECC sonde, 
and from the Dobson spectrophotometer are analyzed by the 
responsible groups, and the observations are provided in 
terms of the various atmospheric parameters. Data selected 
from those reports are entered into the Profile program 
where they are made ready for use with the Rocoz data. 

This section of the document discusses the selection and use 
of the constants and observations. The precise formats for 
the input of data are found in "Profile Program - Input 
Data" and in the chapter on Tape Formats. Within the 
Profile Program, these observations are stored in 
COMMON/ PROFLS/ in arrays which are described in detail in 
the section on Common Areas. 


Constants 


Used in the conversions of units, hypsometric computations, 
estimation of acceleration due to gravity, computation of 
dynamic viscosity, etc., are a number of constants. The 
values for Avogadro's number, the gas constant, the mean 
molecular weights of air and ozone, the volume of an ideal 
gas, Sutherland's constant, and beta are those given in the 
U.S. Standard Atmosphere, 1976. The acceleration of gravity 
at 45 degrees latitude is from the Smithsonian 
Meteorological Tables (R.J.List, ed., 1968, p. 488). 


Avogadro's Number 
gas constant, R* 
acceleration of gravity at 
sea level at 45 deg latitude 
mean molecular weight of air 
mean molecular weight ratio, 
ozone/air 

volume of ideal gas at STP 
Sutherland's constant 
beta 

degrees to radians 


6 . 022169X 10 26 kmol” 1 
8314,32 N m/ ( kmol K) 

9.806160 m sec” 2 
28.9644 kg/kmol 

1.65714 - 

22.4136 nr/kmol 
110.4 K - 

1 . 458x10”° kg/ ( s m K 
0.01745329 rad/deg 
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Dobson Obse-vations 

At the Wallops Flight Facility, total ozone is measured by 
use of a Dobson Spectrophotometer. Occasionally, at remote 
sites, other instruments may be used which observe total 
ozone similar wavelengths, such as the SenTran filter 
photometer or the Brewer Ozone Spectrophotometer. 

To be entered into the Profile Program (Subroutine CARDS) 
are total ozone values observed before and after the Rocoz 
flight; the program refers to these as AM and PM data. The 
date of observation and two values of total ozone are stored 
In words 2-4 of the type -501. pedigree record (PED{x,6)). 
Values of total ozone range from about 200 to 650 Dobson 
Units; a Dobson Unit is intended to be equal to a 
milliatm-cm of ozone but it is recognized that the effects 
of non-ozone absorbers and scatterers may, under some 
circumstances, be inseparable from that of ozone. 

In Subroutine CHECKS the ozone column content as measured by 
the Dobson spectrophotometer is compared to the column 
content indicated by the combined Rocoz-ECC ozone profile. 


Datasonde Observations 

The air temperature profile above about 20 km can be 
obtained from a Datasonde, a small rocket payload consisting 
of a temperature sensor and a telemetry transmitter 
suspended from a decelerator, tracked by radar as it 
transmits temperature data during its descent from the 
vicinity of 60 km. If only one Rocoz flight is planned on a 
given day, one Datasonde is usually scheduled with a launch 
time within one hour of the Rocoz launch. If multiple Rocoz 
flights are scheduled, Datasonde launches are usually 
scheduled shortly prior to the first Rocoz launch and 
shortly after the last Rocoz. Data from the sondes are 
reduced at the Wallops Flight Facility, with scientific 
quality control under the cognizance of F. Schmidlin. The 
reduced data is in terms of air temperature as a function of 
altitude and also includes the 2-sigma error value 
associated with each point. 

To be entered into the Profile Program (Subroutine CARDS) 
are information concerning the flight, the temperature 
profile, and Independently observed air pressure and air 
temperature at a "base" altitude. Identification of the 
flight (an 8-character Flight Number, date and time of 
flight), are stored in words 5-8 of the type -502. pedigree 
record (PED(x,7)). A base altitude and air pressure derived 
from the independent observations made by the ECC sonde are 
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entered and stored in words 9 and 10 of PED(x,7) for use in 
deriving air pressures and densities at the datasonde 
altitudes. The time in minutes (MDEL) between the Datasonde 
and Rocoz launches is entered for use in estimating the 
temperature error. The altitude, air temperature, and 
2-sigma temperature error are entered at integer km 

altitudes, beginning at the highest altitude for which a 
temperature observation is available. The temperatures are 
stored in the remainder of the type -502. pedigree record 

and In as many type -503. records as are necessary. For 
use later within the Profile Program, the base altitude, air 
temperature, and air pressure are stored in DSB in 

C0MM0N/PR0FLS/ , and the altitude, temperature, and 
temperature error are in SONDE, also in C0MM0N/PR0FLS/ . The 
index for the location in SONDE of the lowest altitude level 
of the profile is noted and stored as MXSND. Altitudes are 
stored in km; temperatures in kelvin. 

Because the real interest is in the temperature at the time 
of the Rocoz flight, an error due to atmospheric variability 
is added to the instrumental error. The formula used in 
Subroutine CARDS, from F.Schmidlin, ( J . Geophy s . Res . , 

86,9599-9603 , Oct. 20 , 1981 ) is 

de IT ( Rocoz ) « SQRT ( 0.5 * Terr * Terr ♦ B * MDEL) 

where 

delT(Rocoz) is the one-sigma error in temperature for the 

time of the Rocoz flight, 

Terr is the two-sigma error of the Datasonde 

temperature , 

B is 0.32 above 52.5 km; 0.19 below 52.5 km, 

and 

MDEL is the difference in minutes between the two 

launch times, or a maximum of 100. 

If the flight report does not include estimated errors in 
the temperature, the instrument errors given by Schmidlin 
(see above) may be used. DelT(Rocoz) Is stored in SONDE 
( 5 , x ) and used later in Profile in the estimation of the 
error associated with the expression of the ozone content in 
terms of partial pressure or mixing ratios. 

In Subroutine DATSND, a profile of air pressure and density 
is derived from the datasonde temperature profile and the 
base values of pressure through use of the hypsometric 
equation : 

P(alt) - P(base) * exp (- G * M » (H - HblMOQO) 

(fc * 

where 

P(alt) is the pressure (mbar) at altitude, H 

P(base) is the pressure (mbar) at a base altitude, Hb, 
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G is the acceleration of gravity at the latitude 

of observation and altitude H, 

M is the mean molecular weight of dry air at sea 

level, » 28.9644 kg/kmol, 

R is the gas constant, 

- 8314.32 N m/(kmo1 K) 

T is the average temperature (kelvin) between H 

and Hb, ■ 0.5 (T(alt) ♦ T(base)), 

H is the altitude (km) for P(alt), and 

Hb is the altitude (km) for P(base). 

The acceleration of gravity at the latitude of observation 
is computed according to the Smithsonian Meteorological 
Tables (6th ed., 1 968 , p. 488 , R . J .list ,ed . ) . That is at 
sea level and 45.00 degrees latitude, 

6(45) « 9.806160 m/s*. 

At any other latitude, 

6(lat) - G( 45 ) ( 1 - 0.0026373 cos(2*lat) ♦ 0.0000059 cos*(2*lat)) . 

To adjust for the altitude of observation, 

G(alt) - G(lat) R0*/(RQ*H) 1 , 

where RO is the radius of the earth at the latitude of 
observation and is given by 

RO - 6378.388(1 - 0.0033670 sin* (lat) ♦ 0.0000071 sin*(2*lat)) km, 
from the International Ellipsoid (Basic Physics of the Solar 
System, p.83, V.M. Blanco and S . W. McCuskey , Addison-Wesley, 
1961). 

Densities are then derived from the air temperature and 
pressures as follows: 

D(alt) « (P(alt)MOO . ) * M / (T(alt) * R), where 

D(alt) is the density at each altitude level (kg/m 3 ), 

P(alt) is the pressure at each altitude level (mbar), 

(1 pascal * 100 mbar) 

T(alt) is the temperature at each altitude level 
(kelvin) , 

M and R are the mean molecular weight and the gas constant 
as defined above. 

These data (altitude, temperature, pressure, density, and 
temperature error) are stored in the array, S0N0E (5,70). In 
Subroutine SLANTA, SONDE is used in the computation of slant 
air mass above a given altitude and its rate of change with 
altitude. In Subroutine ALFEFF, SONDE is used in the 
computation of an average temperature, weighted by air 
density for the ten km above a given altitude, for use in 
obtaining the appropriate ozone absorption coefficient. In 
Subroutine 0ZDENS, temperature, temperature error, and 
pressure are stored in the SAVE arrays with the 

corref ponding ozone data and are used in 0ZDENS, C0MPST, and 
RELATD to convert ozone density to ozone partial pressure, 
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the expression of ozone content in terms of mixing ratios 
and layer content, and in the estimation of errors o f these 
derived ozone values. 


ECC Sonde Observations 

One or more ECC sondes are launched near the time of the 
Rocoz rocket launch, usually within an hour, before and 
after. This is a small payload suspended from a weather 
balloon, consisting of a standard rawinsonde with an ECC 
package attached. The data from the electrochemical cell 
(ECC) is inserted into the stream of rawinsonde data, with 
readings initiated by each successive closure in the 
microbaroswitch. The rawinsonde provides data on air 
pressure, air temperature, and humidity, with humidity data 
limited to the lower troposphere. Maximum altitude is 
limited by the burst altitude of the balloon, usually above 
30 km. Altitudes are derived from the observations of 
temperature and pressure via the hypsometric equation; 
recent ECC flights have been radar tracked and these 
altitudes should be more accurate, especially at the higher 
levels. The ECC ozone values, expressed as partial pressure 
of ozone, have their greatest accuracy in the lower 
stratosphere, from the tropopause to no more than about 30 
km. Some profiles show large excursions, usually due to the 
normal "layer-cake" structure of the ozone mixing ratios; 
occasionally there may be an excursion due, apparently, to 
the release of bubbles from the walls of the cell. A 
standard data reduction procedure for ECC's is in place at 
the Wallops Flight Facility, and reports are issued in a 
standard format, providing data at the microbaroswitch 
standard pressure levels as well as at each minute in 
f 1 i ght . 

In Subroutine CARDS of the Profile program, provision is 
made for entering data from two ECC sondes. The flight 
identification (8 characters) and launch date and time are 
stored in type -506. pedigree records. A "tie-on altitude" 
at an integer km level is chosen near the base of the Rocoz 
ozone profile and the corresponding ozone content (atm-cm) 
from ground level to this altitude, readily obtainable from 
the standard ECC data product, are likewise entered into 
Subroutine CARDS and stored in the -506. record. Next are 
entered data for 19 selected standard pressure levels, 
namely at 1000, 850, 700, 500, 400, 300, 250, 200, 150, 100, 
70, 60, 50, 40, 30, 25, 20, 15, and 10 mbars. These data 
(altitude in km, ozone pressure in nanobars, air pressure in 
millibars, and air temperature in kelvin) are stored in the 
array, ECC, starting with the data for the 1000 mbar level 
as well as in the the type -506. and its continuation, type 
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-507. pedigree records. Finally, Subroutine CARDS reads 
data for the ozone column content (atm-cm) from the ground 
up to 250, 125, and 62.5 mbar pressure levels storing these 
three values in ECCLM, for use in Subroutine RELATED. The 
ECC observations are generally the source of information for 
"base" values for use with the Datasonde profile. 

For altitudes below the range covered by the Datasonde 
temperature-pressure profile, the Profile program uses air 
temperatures and air pressures from the ECC data. In 
Subroutine SLANTA air density is used in the computation of 
the slant air mass above various altitude levels. In 
Subroutine ALFEFF, air temperatures and densities are used 
in the computation of a weighted temperature so that a 
temperature -corrected ozone absorption coefficient may be 
obtained. In Subroutine 02DENS, air temperatures and air 
pressures are placed in the type 2xx. data records; the 
one-sigma error in ECC temperature is assumed to be 1.0 
kelvin and the error in ECC pressure is assumed to be 0.7 
mbar. In Subroutine CHECKS, the ozone column content 
between the ground and the "tie-on" altitude is added to 
that inferred from the Rocoz photometer, and the sum 
compared to the column content observed by the Dobson 
spectrophotometer . Finally, in Subroutine RELATD, tables 
are prepared for the user's convenience which include the 
ECC ozone content for standard layers and pressure levels. 


Rocoz Photometer Observations 

Raw data from the Rocoz photometer telemetry tape are 
digitized at the Wallops Flight Facility and processed by 
the RAW, EDIT, MERGE, and SMOOTH programs resulting in a 
SMOOTH tape, which is the input tape to the PROFILE program. 
The radar track of the photometer is likewise digitized at 
the Wallops Flight Facility, processed by the RADAR program, 
and combined with the photometer data in the MERGE program. 
The SMOOTH tape provides photometer current (counts) as a 
function of altitude and related information such as 
calibrations for instrument temperature and battery voltage, 
location of the sun, launch site and time, solar zenith 
angle, and noise content. Details of the tape formats and 
functions of the various programs are given in other 
sections of this document. 

The SMOOTH tape is read by Subroutine TAPE. The Pedigree 
records are stored in the array, PED, for output to the 
Profile tape. The photometer data are placed in the array, 
SMOOTH ( I ,J ,K) , such that all data for S3 have subscript K=l, 
for S2, K-2 , ...for S3/S0, K=4,...and for S1/S0, K=7. In 
Subroutine SECANG, the secant/Chapman-function is computed 
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for the solar zenith angle at each altitude; details of this 
computation are given in the next paragraph. In Subroutine 
OZDENS the photometer values for two-km thick layers are 
interpreted in terms of the ozone density at the center of 
that layer and stored in the array SAVE. The method of 
computing ozone is given in the section on Ozone Algorithm. 
After all the data for the filter under consideration have 
been interpreted in terms of ozone densities, the results in 
SAVE are stored in the corresponding SMOOTH area for use 
later in the Profile program. 

To convert column content from that in the slant path 
between the photometer and the sun to the corresponding 
vertical column content, the slant path content can be 
divided by the secant of the solar zenith angle. This is 
exact for a plane parallel atmosphere, but, for a spherical 
atmosphere, leads to error, particularly for solar zenith 
angles larger than 60 degrees. In Subroutine SECANG, if the 
cosine of the solar zenith angle is less than 0.5, 
approximations to the Chapman function are used for 
computation of the ratio of the vertical path to the slant 
path. The basic equations for the Chapman function and a 
table of values is given by S. Chapman (Proc. Phys. Soc., 
London, 43, 483-501, 1931); this discussion is extended and 
updated tables given in S. Chapman (Proc. Phys. Soc., 

London, B, 66, 710-712, 1953); a simplified form for the 
computation is provided by Fitzmaurice (Appl. Optics, 3, 
640, 1954); the Fitzmaurice approximation was expressed in 
terms of a series expansion by K. Klenk of Systems and 
Applied Sciences Corporation (SASC), and described in their 
Sept. 1980 monthly report on Contract no. NAS5-25346, Task 
no. 41. The SASC expression for the Chapman function is 

Ch = (X it/ 2 ) ^ (a, t ♦ a 9 t 2 ♦ a, t 3 + a. t 4 + a. t 5 ) 

with 1 L J * T 

X = (Re ♦ h)/H , t = 1/(1 ♦ pY), Y = ( X / 2 ) * cos z 

where 

Ch is the value of Chapman function. 

Re is the radius of the earth, km, 
h is the altitude above the spheroid, km, 

H is the scale height of component of interest, km, 

for both ozone and air, set to 5.0, 
z is the zenith angle, 
ir = 3.141593 , 

a, = 0.254829, 

ai = -0.284496, 

a, * 1.421413, 

aj * -1.453152, 

a? = 1.061405, and 

P = 0.327591. 

The above expression is used if the cosine of the zenith 
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angle is less than 0.2 (i.e. for angles greater than 78.5 
degrees). For zenith angles between 60 and 78.5 degrees, a 
more simple approximation is used, namely, 

Ch = (1. - T (I.- 3 T (1 - 5 T) ) ) /cos z 
with , 

T = 0.5 Y*‘. 

In the following table are given values for a precisely 
computed Chapman function at X=1000. for several different 
zenith angles and the values as computed by the the above 
approximations. The precisely computed values are inferred 
from the tables provided by W. Swider Jr. and M. E. Gardner 
(On the Accuracy of Certain Approximations for the Chapman 
Function, AFCRL-67-0468 , August 1967 , 14p, Air Force 
Cambridge Research Laboratories, L.G.Hanscom Field, Bedford, 
Massachusetts ) 



Values of 

Chapman 

Funct i on 


zenith angle 

Chapman 

this 

this 

secant z 

degrees 

X = 100'-. 

o 

o 

o 

tH 

N 

X 

O 

09 

pH 

II 

X 


60. 

i . 9 u 4 

1.992 

1.994 

2.000 

70. 

2.902 

2.899 

2.905 

2.924 

80. 

5.590 

5.598 

5.643 

5.759 

90. 

39.648 

39.633 

44.874 



Note that X = 1 4 8 0 . corresponds to an altitude of 40 km and a 
scale height of 5 km. 


Radar Tracks 

Data for the radar track are digitized at the Wallops Flight 
Facility, smoothed by the Code 970 Eclipse computer (MESUP 
tape), converted to metric units by the RADAR Program, and 
merged with the Rocoz photometer data in the MERGE Program. 
L ror in altitude for data from the WFF radars is assumed to 
be 10 m and is used in Subroutine OZDENS in the estimation 
of error in d ( 1 n I ) / d h . Altitude error for observations at 
remote sites is probably larger . 


Model Ozone Profile 

A model ozone profile is provided by Subroutine MODEL for 
use in Subroutine SLANTA in its estimate of the ozone column 
above the top of the ozone profile from the Rocoz data under 
consideration, and, in Subroutine CHECKS, for comparison 
with the observed Rocoz ozone profile. 
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Number densities (molecules m are from the mid -latitude , 
Northern hemisphere model provided in the U.S. Standard 
Atmosphere, 1976, based on an analysis by A. Krueger and R. 
Minzner and are stored in a DATA statement. Values are 
interpolated so as to provide an ozone number density at 
each km from 11 to 70 km. Subroutine MODEL assumes an 
overburden above 70 km of 5.35x10 atm-cm and infers the 
overburden at lower altitudes from the model density 
profile. To obtain the column content, the density is 
multiplied by 1.0x10 x 22.4136 liters per mole / 
6.022169x10 (Avogadro's number) and added to the column 
content at the previous altitude level. The total column 
above 0 km is 0.345 atm-cm. A more elegant method would be 
to use tables derived from satellite data to obtain a "first 
guess" profile appropriate for the latitude, season, and the 
column content as observed by the Dobson spectrophotometer . 
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Effective Ozone Absorption Coefficient 


Computations 

Given values for the content of ozone and air in the slant 
path to the sun, characteristics of the filter, and the 
necessary reference tables. Subroutine ALFEFF returns values 
for the effective ozone absorption coefficient and the 
effective scattering coefficient. 

The filter shape may be expressed in terms of its 
transmission as a function of wavelength. Such measurements 
have been made, in the past, at the Naval Weapons Center 
(NWC) at China Lake, California, and are now done at the 
Wallops Flight Facility. These tables may be read by 
Subroutine SHAPE and used in Subroutine ALFEFF to compute 
Aeff. The formulae are similar to those used by NWC as 
given by h. E. Hills and C. A. Flanagan (Calculation of 
Effective Ozone Absorption Coefficients for a Rocketborne 
Ozonesonde, NWC TP 5904, 41 pp. Naval Weapons Center, China 
Lake, Calif., March 1977). That is, 


Aeff = J. In sum 1 0 ( L ) Q(L) S(L) F ( L ) exp(-A(L) X -B ( L ) M) dL 
X sum 1 0 ( L) Q(L) S(L) F ( L) exp(-B(L) M) dL 


or for slant path ozone less than l.xlO" 5 atm-cm and slant 
air mass less than 0.01 atm, i.e., no overburden. 


Aeff = sum I0(L) Q( L) S(L) F ( L) A ( L ) dL 
sum 1 0 ( L) Q(L) §( l5 F|lJ dL 

where 

L is wavelength, at 0.1 nm intervals, 

10 is the flux of the sun at zero optical depth, 

Q is the detector respons i vity , 

S is the diffuser transmission, 

F is the filter shape, 

A is the ozone absorption coefficient, per cm, base e, 
X is the slant path ozone overburden, atm-cm, 

B is the scattering coefficient, per atm., base e, and 
M is the slant path air mass, atm. 


The ozone absorption coefficient is a function of 
temperature, particularly for wavelengths longer than about 
270 nm. In Subroutine ALFEFF, a weighted temperature is 
computed, dependent on the air density for ten km above the 
altitude under consideration, and the ozone absorption 
coefficient interpolated at each wavelength from tables of 
the ozone absorption coefficient at several different 
temperatures . 
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To obtain an effective scattering coefficient, a similar 
equation is used, that is, 

Beff = -1 In sum I0(L) F ( L) Q(L1 S(L) exp(-B(L) M) dL 

“ff sum 16) L) F(L) Q(L) I(L) dL — — . 

The Solar spectral flux, Rayleigh scattering coefficients,' 
ozone absorption coefficients, diffuser transmission, 

detector responsi vity, and filter transmission tables are 
stored on disk and read by Subroutine INARRY. The diffuser 
transmission and detector responsivity are assumed to be the 
same for all photometers. The filter transmission is 
measured for each of the four filters in each photometer, 

and must be supplied for each flight. The filter shapes, 

for the flights up through Rocoz no. 315, are expressed only 
in terms of coefficients, and are read by Subroutine CARDS. 
For later flights, detailed filter shapes are also 
available, and may be read by Subroutine SHAPE. 

Subroutine INARRY expects the disc data sets to have an 
initial 80-byte header followed by 1000 real numbers 
corresponding to wavelengths of 240.1 to 340.0 nm at 0.1 nm 
intervals. The header consists of 64 characters followed by 
4 real ( R* 4 ) words. The real numbers are the range of 
wavelengths and the temperature, as applicable. The first 
part of the character string should identify the contents, 
including the parameter, the observer, and the date, as 
appropriate . 

If McBride coefficients are available (not equal to zero), 
then Subroutine ALFEFF uses an approximation developed with 
NWC describing the dependence of Aeff on the ozone content 
of the slant path, namely, 

Aeff - A0 ♦ 2. A1 X ♦ 3. A2 X 1 
where 

X is the slant path ozone overburden, atm-cm, 

A0 is the effective ozone absorption coefficient for no 
overburden, and 

A1 and A2, as well as A0, are coefficients describing 
the change as a function of overburden. 

The three coefficients, A0, Al, and A2, are obtained by 
computing Aeff for various overburdens and fitting the 
results to a polynomial, assuming that the effects of 
scattering on Aeff are negligible above about 20 km. The 
coefficients are part of the card (DATA5) input and are 
computed prior to running the PROFILE program. 
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The effective scattering coefficient, Beff, is set equal to 
value of B provided by the NWC analysis and is assumed to be 
constant over the range of air masses encountered in the 
reduction of the Rocoz observations. It is estimated that 
that the error due to this assumption is no more than 2% for 
altitudes above 20 km. The NWC analysis also provides the 
center wavelength and width at half height for each filter. 
These may be used, with the assumption of a gaussian filter 
shape, to estimate the effect of the temperature dependence 
on the Aeff. 


Solar Flux 

For both the solar flux and the ozone absorption 
coefficient, several different sets of data are available. 
For flights prior to no. 315, a "McBride" solar flux was 
used at the Naval Weapons Center (China Lake, California) in 
the computation of effective absorption coefficients for 
each filter. This is based on the average values for 0.1 nm 
intervals of the Furukawa compilation (Furukawa, P.M., 
P.L.Haagenson, and M . J . Scharberg , A composite, High- 
Resolution Solar Spectrum from 2080 to 3600 A, NCAR-TN-26, 
55p, 1967, National Center for Atmospheric Research, 
Boulder, Colorado) and is available as a disc data set for 
use with the Profile Program. 


Rayleigh Scattering Coefficient 

The Rayleigh scattering coefficient as a function of 
wavelength was likewise assembled by Dr. McBride and is 
available as a disc data set for use with the Profile 
Program. 


Ozone Absorption Coefficients 

The effective absorption coefficients computed by Dr. 
McBride were based on ozone absorption coefficients measured 
by E.C.Y.Inn and Y. Tanaka (Ozone Absorption Coefficients in 
the Visible and Ultraviolet Regions, Ozone Chemistry and 
Technology, Advances in Chemistry Series No. 21, pp263-68, 
1959) for wavelengths shorter than 272 nm and from 
measurements at -44 C. reported by E. Vigroux (Contribution 
6 l'6tude experimental de 1'absorption de l'ozone, 
Ann.Phys., 8, 709-62, 1953) for wavelengths longer than 290 
nm, with an appropriate interpolation at intermediate 
wavelengths. This data set is available on disc for use by 
Subroutine INARRY. INARRY expects three sets of ozone 
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absorption coefficients at three different temperatures 
(such as are expected to be available from the work of 
A. Bass). In Subroutine ALFEFF interpolation is made to 
obtain the ozone absorption coefficient corresponding to the 
air temperature weighted according to air density. 


Diffuser Transmission 

The transmission of the three stacked diffuser plates as a 
function of wavelength was measured a number of years ago, 
and the resulting curve has been used for all filters 
characteri zed by Dr. McBride. This is available as a disc 
data set for input to the Subroutine INARRY. 


Detector Responsivity 

The response of the detector was obtained many years ago 
(from the manufacturer ?), was used by Dr. McBride in his 
analyses, and is available as a disc data set for input to 
the Subroutine INARRY. 


Filter Transmission 

The shape of each filter (with the nickel sulfate 
hexahydrate crystal) is measured by means of a Cary 14 Dual 
Beam Spectrophotometer, and the measured transmission is 
used by Subroutine ALFEFF in the computation of the 
effective absorption coefficient. 

For the filters calibrated in Dr. McBride's laboratory at 
the Naval Weapons Center(NWC), an NWC analysis provides the 
center wavelength and half-width, the effective absorption 
coefficient with no overburden (AC), and the coefficients, 
A1 , A2 , and B . 

For the filters calibrated at the Rocoz Test and Calibration 
Laboratory at the Wallops Flight Facility, Subroutine SHAPE 
reads the data set containing the measured spectral 
transmission, and Subroutine ALFEFF computes effective ozone 
and scattering coefficients as a function of air temperature 
and air and ozone overburdens. 
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Ozone Profile 


Basis of Ozone Algorithm 

The technique for converting the observed solar flux as a 
function of altitude to ozone content as a function of 
altitude was developed by Arlin Krueger (NASA, Goddard Space 
Flight Center) for use with his Rocoz photometer and was 
described in detail by Eugene Shaffer (Systems and Applied 
Sciences Corporation) in a draft report. This discussion is 
based primarily on that report, but adds the details of its 
implementation in the Profile Program. 

The reduced Rocoz data consist of a flux of ultraviolet 
light in four different wavelength bands as a function of 
altitude. The wavelength bands are chosen over a range of 
ozone absorption coefficients. To compute the ozone 
profile, the following assumptions are made: 


1. Dust particles cause a negligible attenuation of 
the direct sunlight. 

2. No other absorbers are present in the stratosphere 
for the four wavelength bands. 

3. Above 20 km, the diffuse component of the sunlight 
can be ignored. 


With the above assumptions in mind, the intensity of the 
sunlight as a function of altitude can be written, using 
Beer's law, as 


I(L,h) = I { L ,• ) exp [ A( L) u ( h ) - B(L)m(h)] 


( 1 ) 


Where: I(L,h) is the intensity at height h, 

I ( L , • ) is the intensity at top of the atmosphere, 

A ( L) is the ozone absorption coefficient, 
u(h) slant path ozone from h to the sun, 

B(h) Rayleigh scattering coefficient, 
m(h) optical air mass between h and the sun, 

L is wavelength, and 

h is altitude. 

Equation (1) can be rewritten as 


In I ( L ,h) - In I(L,«) * -A(t)u(h) - B(L)m(h). 

An expression for the ozone concentration at each height 
Interval can be derived by differentiating the above 

expression with respect to altitude. 

d In I(L.h) - -A ( h ) du(h) - B(h) dm ( h ) (2) 

dh dh dh 


11-47 



PROFILE PROGRAM 


The monochromatic ozone absorption coefficient and Rayleigh 
scattering coefficient must be replaced by effective 
coefficients which are now a function of the slant path 
total ozone and the air mass amounts. 

A ( L ) is replaced by Aeff(u,m) 

B(L) is replaced by Beff(u,m) 

Effective coefficients are used because the optical filters 
have bandwidths wide enough that the wavelength dependence 
on various input parameters must be recognized. The slant 
path ozone amount, (u), and the optical air mass amount, 
(m), modify the spectral shape of the solar flux as a 
function of height. 

The effective coefficients are substituted into equation (2) 
to yield 


d In I(L,h) - -Aeff du - d(Aeff) du u - Beff dm (3) 
'dh du iJh 3K 


Equation (3) is rewritten to represent du/dh, the total 
amount of ozone in the height interval. 


du 

3K 


! L _Qn H LJiIl 

Tell ♦ d (Ie77) u { (fn 
du 


♦ Beff dm) (4) 
Hn } 


Finite differences may be substituted for the derivatives, 
and equation (4) now can be solved iteratively for ozone 
amounts at any height level in the range of the Rocoz data. 

The implementation of this solution in the Profile Program 
makes use of the following: 

- The magnitude of the solar flux at the wavelengths of 
interest at any given altitude is taken to be 

proportional to the current from the photodiode 
detector as adjusted by the Rocoz electronics for 
variations in the flux at longer wavelengths entering 
the photometer, that is, I is proportional to the 
number of counts transmitted by the telemetry system 
for any particular filter. The constant of 

proportionality is not important since the difference 
of the logarithms is the parameter that is used. 
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V 


One way to represent the effect on the ozone 
absorption coefficient is to use a quadratic function: 

Aeff(u,m) » AO ♦ A1 u + A2 u* * B1 m ♦ B2 m* 

At typical rocket heights (20-60 km) the fourth and 
filth terms in equation (3) are negligible due to the 
small air mass and are Ignored in this discussion. 
The effective Rayleigh scattering term can be 
considered a linear function of the optical air mass. 
For heights of 20 to 60 km we represent the effective 
Rayleigh scattering as a constant, which for the 
filters in use here will have an error of no more than 
2X. The constants in the quadratic equation may be 
calculated through use of the measured filter shape, 
solar spectrum, ozone absorption coefficients, 
Rayleigh scattering coefficients, and the spectral 
response of the Rocoz photometer, together with 
various assumed ozone and air overburdens. These 
constants have been calculated by Dr. William R. 
McBride for the filters that have been calibrated at 
the Naval Weapons Center at China Lake, California, 
and are available for use ir. the analysis of the Rocoz 
data for flights through 1982. 

It is noted that dm/dh is simply (air density x secant 
solar zenith angle). 

The Profile Program will process the Rocoz photometer 
data in terms of the observed count for each of the 
four filters (S3, $2, SI, and SO) or in terms of the 
ratio to the SO filter (S3/S0, S2/S0, S1/S0). The use 
of ratios is preferable if there is appreciable common 
mode noise and/or if there is a change in the 
responsivity of the photometer as a function of time 
in in f 1 i ght . 

The range of altitude for each filter is chosen prior 
to running the Profile Program; it generally 
corresponds to optical depths between 0.05 and 3.0 and 
depends on the availability and quality of data, that 
is, both radar and telemetry data must be available, 
telemetry dropouts must be minimal, the attitude of 
the Rocoz must be reasonably stable, and the solar 
zenith angle should be less than 60 degrees. Data for 
solar zenith angles up to about 80 degrees can be 
processed, but the accuracy of the resulting ozone 
densities becomes increasingly uncertain. 
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Implementation of the Ozone Algorithm 


The Rocoz data is interpreted in terms of ozone densities by 
the Profile Program in subroutine OZOENS. The evaluation 
involves three altitude levels, spanning a total of 2 km, 
namely, top (T), center (C), and bottom (B), with the 
differentials evaluated for a layer centered at C. In the 
initial evaluation, estimates are used for ozone densities 
and ozone overburdens. After the initial evaluation of the 
ozone density at C, these estimates are updated and the 
ozone density at C recomputed until the change in the ozone 
density compared to its previous value is less than 0.1 X. 

The steps necessary to interpret the input data are as 

fol lows : 

1. The top (ALTT) , center ( ALTC } , and bottom ( ATLTB ) 
altitudes are identified, beginning with ALTT in km as 
provided in the array TOP, and ALTC and ALTB one and two 
km lower, respectively. After the determination of the 
ozone density at ALTC, ALTT, ALTC, and ALTB are “*c., 
decremented by one km, ending with the layer for \ 'c h 
ALTB corresponds to the value given in the array BAS... 

2. SLANTA is called to provide the slant air mass (ASLTT, 

ASLTC, ASLTB) above each of the three altitudes and its 

rate of change evaluated at the center altitude (DMDHC). 
Slant air mass is in atmospheres. DMDHDC is in 

atmospheres per km. 

3. The rate of change of the natural logarithm of the 

intensity, d (ln(I))/dh, (DLIDH) is evaluated from the 
observed intensity (or ratio to SO) corresponding to 
ALTT and ALTB, the top and bottom of the three layers 
under consideration. 

4. Estimates of the slant ozone overburden (UT, UC, U B ) are 

made for each of the three altitudes. Initial estimates 
are taken from the from the model ozone profile 
contained in the array STD. Ozone values Inferred from 
the rocket flight are used In the estimates for 

succeeding levels. Ozone overburden is in atm-cm. 

5. Values at each altitude for the effective ozone 
absorption coefficient (ALFAT, ALFAC, ALFAB) and the 
effective Rayleigh scattering coefficient (BETAT, BETAC, 
BETAB) are obtained from subroutine ALFEFF. 
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6. The slant ozone content between ALTT and A..TB In atm-cm 
(DELU) Is computed using equation (5) < nd saved as 
(PDELU) . 

7. DELU is used to recompute the slant ozone for the bottom 
layer (UB) t ALFAB and BETAB are recomputed, ALFAC is 
taken as the average of ALFAT and ALFAB, and BETAC is 
the average of BETAT and BETAB. 

8. Steps 6 and 7 are repe^ed until DELU is within 0.1% of 
PDELU, usually requiring but one iteration. 

9. Ozone density at ALTC is then 0.5 * DELU /(secant of the 
solar zenith angle), in atm-cm/km. 

10. The ozone density and a number of related values are 
saved in the type lxx. and 2xx. Profile data records. 

11. Steps 1-10 are repeated at the next lower altitude level 
until the range for a particular filter is completed. 

12. Steps 1-11 are repeated for the next filter until all 
seven possibilities have been considered. 


Profiles from Each Filter 

As the ozone density is obtained at each altitude level from 
each filter (or from its ratio to the SO data) two data 
records are created and are included on the Profile Tape. 
(See the type lxx. and 2xx. data records described in the 
chapter on Tape Formats.) The first record contains 
relevant parameters taken from the Smooth Tape and various 
elements calculated in the course of obtaining the ozone 
density. The second record is devoted to atmospheric 
parameters, that is, ozone values, air pressure and 
temperature, and associated errors. Selected values from 
these records are printed as tables. 
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Composite Rocoz Profile 

The Rocoz photometer contains five filters, referred to as 
SO, SI, S2, S3, and Compensation, with S3 having the 
shortest center wavelength, S2, SI, and SO at longer 
wavelengths, and the Compensation filter at an even longer 
wavelength for which the ozone absorption coefficient is 
very small. The signal from the detector with the 
Compensation filter is used as the reference in an automatic 
gain control system to continually change the gain of the 
amplifier for the detector associated with the S0-S3 
filters, compensating for changes in total light input such 
as those due to changes in the direction of the axis of the 
photometer with respect to the sun. 

Ozone profiles are derived independently from the data from 
each of the four filters (SO, SI, S2, S3). Each profile 
covers an altitude region of 10 to 25 km, usually with some 
overlap in altitude coverage. In Subroutine COMPST these 
independent profiles are combined to form a single profile 
attributable to this Rocoz flight. To obtain this composite 
profile, all available data for a particular altitude level 
are weighted by the inverse of the associated error and 
averaged. Ozone densities derived from filter ratios are not 
included . 

For the convenience of the user, the ozone ..density profile 
is expressed in both atm-cm/km and number/nr and is collated 
with air temperatures and pressures based on the Datasonde 
and the ECC sonde. In Subroutine CHECKS, the column content 
(atm-cm) at each altitude level is computed, based on an 
initial column content above the top level of the ozone 
profile taken from the ozone model (Subroutine MODEL), and 
integrating downward. 

In Subroutine CHECKS the Composite data are used in a series 
of consistency checks before filling the remainder of the 
array with computed partial pressures and mixing ratios. 
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Consistency Checks 

To aid in estimating the quality of the data from a Rocoz 
flight and related observations, the data are checked for 
internal consistency. These checks include (1) comparison 
of the profiles from the various filters with each other, 
(2) comparison of the composite profile with a model, (3) 
comparison of the ozone column content measured by the 
Dobson spectrophotometer with the sum of the columns 
observed by the Rocoz photometer and the ECC sonde, (4) 
comparison of the slope deduced in the Smooth Program over a 
variable height interval with that deduced in the Profile 
Program over a 2 km interval, (5) comparison of the 
variations of the Compensation channel with those expected 
due to changes in solar zenith angle and overhead air mass, 
and (6) the drag coefficient as a function of altitude and 
Reynolds number. 

The following discussion gives the basis for each check and 
the significance of the results. 


Overlap Analysis 

One way to examine the accuracy and precision of the ozone 
densities is to compare the profiles obtained from the 
different filters. In Subroutine COMPST data for the various 
filters are transferred to a single table; Subroutine CHECKS 
calls Subroutine OVERLP to systematically compare the data 
from the various filters and to compare data from a 
particular filter to data derived from the ratio of that 
filter to the SO filter. For each of the 21 possible 
combinations, the two altitude profiles are compared by 
fitting the pairs of points to a straight line (Subroutine 
UNFIT), and computing its slope { B ) and zero intercept (A), 
i . e. , 

03(filter y) = A + 6 03(filter x). 

One would look for a zero intercept that is small compared 
to the ozone value at the maximum altitude of the set of 
points, for a slope that is near 1.0, and for a correlation 
coefficient (R) that is near 1.0. An examination of the 
table may reveal whether one of the filters has 
substantially changed its spectral transmission since its 
calibration. Values of B that are markedly different from 
1.0, but similar from flight to flight may indicate that one 
or more of the spectral data sets (solar flux, absorption 
coefficients, etc.) is not correct. One would expect that 
the values for A, B, and R for a filter paired with its 
ratio (e.g. S3 compared to S3/S0) would be very close to 
ideal . 
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Comparison with Model 

The ozone density and overburden (vertical column) as a 
function of altitude from the Rocoz Composite profile are 
compared with the values provided by Subroutine MODEL. 
Currently the model is for mid-latitudes and is from the 
U.S. Standard Atmosphere, 1976. It would desirable to 
upgrade Subroutine MODEL to provide a profile that is 
typical of stated latitude and season, such as might be 
might be obtained from the analysis of global ozone data 
sets. Subroutine CHECKS computes the differences and 
prepares a table of the results. 

One would expect most data from Wallops Island to be 
relatively close to the mid-latitude model, within the 
variability stated in the tables given in the U.S. Standard 
Atmosphere, 1976. 


Total Ozone 

The total ozone (vertical column content) is observed by the 
Dobson Spectrophotometer and may also be deduced from a 
combination of the ECC and Rocoz Composite profile. The ECC 
values are generally considered to be relatively accurate 
below the ozone density maximum. The Rocoz profiles are 
better at the higher altitudes. 

From a preliminary examination of the data, a "tie-on" 
altitude is chosen, usually between 18 and 25 km; the 
column content from the ground to the "tie-on" altitude is 
obtained from the ECC data report and entered into the 
Profile Program (Subroutine CARDS). In Subroutine CHECKS, 
the column content above the "tie-on" altitude is obtained 
from the Rocoz Composite ozone profile. Errors are assumed 
to be 7% for the Dobson total ozone; 10% for the ECC 
segment, and as derived in the Profile Program for the Rocoz 
segment. The two segments are added, the associated errors 
derived, and the difference from the Dobson total is noted. 

The two observations of total ozone should agree within 
their estimated errors. If the difference is larger, it may 
indicate temporal and spatial variability of ozone 
(particularly if maps from TOMS, the total ozone mapping 
spectrometer on Nimbus-7, indicate a substantial gradient in 
that area). A large difference may also indicate that one 
or more of the three systems involved did not perform as 
expected and has as yet unidentified errors. This check is 
sensitive only at altitudes where the ozone density is 
large, that is, in the lower stratosphere . 
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Estimations of Slope 

In the Smooth Program, all points in the vicinity of each 
integer altitude are considered, and a value of the 
intensity (I, number of counts) representative of that level 
is derived. The number of points considered range between 
100 and 800; the altitude interval ranges between 1 and 
9 km. In the process, a line of the form, ln(I)= A + 
B(h-hb), is used, and the intercept (A) at the base altitude 
(hb) and the slope (B) are calculated; A and B are then 
used to obtain the intensity at the altitude level under 
consideration . 

In Subroutine 0Z0ENS of the Profile Program, d ln(I)/dh is 
obtained by taking differences over a 2-km interval. For 
ideal data, the two values should be the same, at least 
within the error associated with each. Larger differences 
are probably related to the differing altitude resolution, 
but should be examined individually. 


Variation of the Compensation Channel 


A fraction of the light entering the Rocoz instrument is 
diverted through a Compensation filter to a separate 
detector. This is a relatively broad-band filter centered 
at about 375 nm, where the ozone absorption coefficient is 
several orders of magnitude less than that for the SO 
filter. Variation of the signal from the detector is due 
primarily to changes in the angle of the sun with respect to 
the optical axis of the photometer, but is also affected by 
the amount of light scattered out of the optical path by 
molecules and aerosols. In the Smooth Program, the value of 
the Compensation signal is averaged over 99 filter 
rotations. In Subroutine CHECKS of the Profile Program, the 
average value (after taking into account changes in solar 
zenith angle) between 35 and 50 km is determined; the 
expected value at lower altitudes is computed, again taking 
into account the changes in solar zenith angle and, in 
addition, the loss of light expected due to Rayleigh 
scattering. The fractional loss is estimated by the method 
developed by R.W. Fraser (NASA/GSFC, quoted (incorrectly) in 
Thomas, R.W.L., W. A. Pearce, A. H. Holland, and D.U. Wright, 
Appl. Optics, 21, 2436-2441, 1982). The correct formula is: 


fc = 

where 

fc is the 
z is the 
h is the 
Differences 


0.376 (sec z) 0 * 874 io" h/16 


fraction of light scattered out of the path, 
solar zenith angle, and 
altitude in km. 

between the predicted and the observed values of 


r 
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the average Compensation signal may be due to scattering by 
aerosols, changes in temperature of the photometer, or large 
(more than 60 degrees) solar zenith angles. 


Drag Coefficient 

The photometer is suspended from a Starute decelerator, 
deployed after ejection from the rocket body. The drag 
coefficient of the Starute and photometer can be computed 
using 

C * (g - a) / (2 d v* A) 

where 

C is the drag coefficient, 
g is the acceleration due to gravity, 
a is the observed vertical acceleration, 
d is air density, 

v is the observed vertical velocity, and 
A is the effective area of the Starute. 

Because drag is a function of Reynolds number, this also is 
calculated and listed. That is, 

Re = v D d / m 

where , 

Re is the Reynolds number, 

D is the character i st i c size, 
m is the dynamic viscosity, and 
d and v are as defined above. 

The dynamic viscosity is computed by use of the formula 
given in the U.S. Standard Atmosphere, 1976, (p.19), 

m = B T 372 / (T + S) , 

where, _ 6 i 

B is a constant, equal to 1.458 x 10 k g / ( s m K 1 ) , 

T is temperature in kelvin, and 
S is Sutherland's constant, equal to 110.4 K. 

The drag coefficient should vary smoothly as a function of 
Reynolds number. The observed values increase from about 1 
near Apogee to about 1.4 near 20 km. This increase is 
similar to that for a sphere at similar Reynolds numbers. A 
substantial deviation from the smooth increase with 
decreasing Reynolds numbers could indicated problems in the 
altitude versus time data; A decrease in the magnitude of 
the drag coefficient could be due to a damaged Starute. 
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Deri ved Ozone Values 


The Rocoz data as analyzed in this set of programs is 
interpreted in terms of ozone density as a function of 
geometric altitude. Although expressed at 1 km intervals, 
the smoothing interval varies from 3 km at low altitudes to 
about 9 km near the top of the profile. The density profile 
is integrated from the top down to obtain the column content 
as a function of altitude. 

For the convenience of the user, the ozone content is also 
expressed in other commonly used units, that is, as partial 
pressure, as mixing ratio, as layer content, and with 
altitude expressed in terms of pressure level. The 
associated error increases because it is necessary to use 
the observed air temperatures and pressures in the 
transformation of units, and the errors in such observations 
must be considered. 


Column Content 

The slant ozone (ozone in the path between the photometer 
and the sun) and the ozone overburden (vertical column above 
the photometer) are computed in Subroutine OZDENS and given 
to Subroutine ALFEFF for use in the computation of the 
effective ozone absorption coefficient. The initial value 
(at the highest altitude for which there is useful 
photometer data, that is, 1 km above the top altitude of the 
Rocoz ozone profile) is obtained from the model profile 
provided by Subroutine MODEL. Values of slant ozone at 
lower levels are obtained by adding the observed layer 
content. The vertical column for each level is obtained by 
dividing by the secant/Chapman-function of the solar zenith 
angle at that level, as provided by Subroutine SECAN6. 
These values are computed independently for each filter, 
i.e., the initial slant ozone used in deriving an ozone 
profile from the S2 filter is taken from the model profile, 
not from the profile obtained by use of the S3 filter. The 
values used in the computations are included in the type 
lxx. and 2xx. data records created by the Profile Program. 

After ozone profiles are obtained corresponding to each of 
the filters of the photometer, a Composite Rocoz ozone 
profile is generated by Subroutines COHPST and CHECKS, as 
described in an earlier section of this document, and 
included in the type 300. data records. The initial value 
of the ozone overburden here Is that from the ozone model, 
multiplied by the ratio of the ozone density in the 
Composite profile to that in the model profile. The density 
for each lower level is obtained by assuming the density in 
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the 1 km layer is the average of the density at its top and 
base, and adding the layer content to the overburden already 
available for the top of the layer. These overburdens are 
the ones used by Subroutine CHECKS for comparison with the 
model overburdens and the Dobson total ozone observations. 


Partial Pressure and Mixing Ratios 


The ozone algorithm as implemented in Subroutine OZDENS 
expresses the ozone density at a given altitude level in 
terms of atm-cm/km. The following formulae are used in 
Subroutines OZDENS, CHECKS, and RELATD to express the ozone 
density in other units: 

0(f) * 0(ack) 1 x 1 0 ” 5 A / V 

0 ( n b ) - 0(ack) 0.1 R* T /V, 

O(vmr) = 0 ( n b ) / P, 

O(mmr) = 1.65714 O(vmr), 
where » 

0(f) is ozone density in number/m , 

O(ack) is ozone density in atm-cm/km, 

0 ( n b ) is ozone density in nanobars, 

O(vmr) is ozone volume mixing ratio, parts per million, 
O(mmr) is ozone mass mixing ratio, parts per million, 

P is air pressure, millibars, 

A is Avogadro's number, 6.022169x10^° molecules/kmol , 

R* is universal gas constant, 8314.32 N m / ( K kmol), 

T is air temperature, kelvin, and - 

V is volume of ideal gas at STP, 22.4136 nr 5 / kmol. 


Content of Standard Pressure Layers 

To facilitate the comparison of the Rocoz/ECC profiles with 
those from other sources. Subroutine RELATD prepares a table 
of layer contents, based on the Composite Rocoz profile and, 
at the lower altitudes, on the concurrent ECC sonde. Layer 
boundaries are expressed in mi 11 i atmospheres beginning at 
1000 matm with the upper boundary at one-half the pressure 
of the lower boundary. The topmost layer is from 0.24 - 0 
milliatm. The layer's ozone content (in mil 1 iatmosphere-cm) 
is the difference of the ozone overburdens at the layer 
boundaries, as interpolated from the data in the Rocoz 
Composite ozone table. The ozone content for layers for 
which the Rocoz profile is not complete is set to 0. 

At altitudes in the troposphere and lower stratosphere, 
layer content is derived from the ECC sonde that was 
launched near the time of the Rocoz launch, using the column 
contents given in the ECC tabular listings. 
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Ozone Values at Standard Pressure Levels 

At 33 standard pressure levels ranging from 1000 mbar to 
0.05 mbar. Subroutine RELATD displays ozone number densities 
interpolated from the Rocoz Composite profile and the ECC 
sonde, air temperatures from the Datasonde and the ECC 
sonde, geometric altitude at each pressure level, and the 
ozone mass mixing ratio. Where the various data sets 
overlap, two sets of temperature, ozone density, and mass 
mixing ratio are included. The standard pressure levels are 
at decimal multiples of 10, 15, 20, 30, 40, 50, and 70, and 
also include 250 and 650 mbar. 
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Estimation of Errors 


There are many sources of uncertainty which affect the ozone 
profiles obtained through use of observations from the Rocoz 
photometer, and the associated Datasonde, ECC sonde, and 
Dobson spectrophotometer. Besides errors related to the 
Rocoz photometer, there are temperature sensor errors, 
altitude (radar) errors, pressure (microbaroswi tch or 
hypsometer) errors, and uncertainties due to atmospheric 
variability (air temperature, air pressure, and Rocoz 
photometer observations are not from the same vehicle), 
variability in aerosol content, and solar variability. 

The error in the ozone density derived from the Rocoz 
photometer depends on the accuracy of the various tables 
used in the computation of the effective ozone absorption 
coefficient: the ozone absorption coefficient and its 
temperature dependence, the solar flux, the Rayleigh 
scattering coefficient, and the measured filter 
transmission, diffuser transmission, and detector 
respons i vi ty . The error in the ozone densities also depend 
on the effectiveness of the compensation scheme, detector 
noise, noise introduced by the in-flight digitization, noise 
introduced by the telemetry and recording systems, the 
accuracy of the radar tracks, and approximations and 
assumptions made in the course of the data reduction. 

The Profile Program attempts to evaluate many of these 
errors. However, it does not characterize errors in the 
computation of the effective ozone absorption coefficient, 
except as can be done by comparing the ozone densities from 
the different filters. Most errors are expressed as percent 
and all represent the one-sigma value. 


Irradiance as a Function of Time 

In this scheme of data analysis, the observed irradiance is 
always expressed in arbitrary units; no attempt is made to 
convert the arbitrary units to absolute fluxes. Ozone 
content is obtained by using the ratio of the radiance 
observed at the top of a chosen layer to that observed at 
the base of that layer and hence is independent of the 
absolute values of the radiances. The unit of irradiance is 
the number of counts (0 - 1000) at the input of the Rocoz 
telemetry transmitter and represents the solar flux received 
on a horizontal plane located at the altitude of the 
photometer. 
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The errors in irradiance are handled in three different 
ways, depending on the time scale. The first group of errors 
have periods of well under one second, the second group have 
periods on the order of 1 to 100 seconds, and the last group 
are drifts over the 20 minutes of the flight. 

The first group of errors relate to the determination of the 
value of a single observation. The output of the photodiode 
is sampled at precisely the time at which the filter is 
fully in the optical path and is held electronically for one 
rotation of the filter wheel. During approximately 600 ms 
required for one rotation, the signal on hold is sampled 
about 55 times. The design of the electronics restricts the 
decay of the signal to less than one count; the 
analog-to-digital conversion is designed to have an 
accuracy of + 1 count. The extent of the decay and the 
short-term random component of the a/d conversion can be 
checked by a visual inspection of the listing of the 
telemetry data. Noise in the photometer electronics may 
correspond to several counts and noise introduced during the 
telemetry transmission and recording may be as much as full 
scale. In Subroutine AVG of the Edit Program, it is 
required that at least half of the samples be within 6 
counts of each other, the remainder are discarded, and an 
average is taken of the selected values. Otherwise a dummy 
value (-99.) is used to indicate that no value is available 
for this filter for this rotation. No estimate is made of 
the error attributable to these sources alone. 

The next group of errors are of somewhat lcwer frequency and 
are due primarily to low frequency components in the noise 
inherent in the photodiode and its current amplifier, to 
noise in the compensation circuitry, and possibly due to 
less than perfect optics such that azimuth dependent 
components may be present in the observed irradiance. In 
the Smooth Program, a value of the irradiance at each 
integer altitude is obtained from consideration of a number 
of observations in the vicinity of that altitude. The 
minimum number of points assembled by Subroutine SELECT is 
100; the minimum altitude range is 1 km; at high altitudes 
the 100 points correspond to an altitude interval of about 8 
km; at low altitudes several hundred points may be needed to 
cover one km altitude. In Subroutine LOGFIT, a straight 
line of the form. In I * A ♦ B (H - Hb) , is fitted to this 
array of points and values of A and B (with standard 
deviations) are obtained. Here, I is the irradiance in 
terms of counts, H is the altitude of the point under 
consideration, and Hb is the lowest altitude in the interval 
under consideration. The irradiance at the integer altitude 
level is then computed using the computed values of A and B. 
This value of irradiance, A, B, the altitude bounds, and the 
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standard deviations associated with A and B (sigmaa and 
sigmab) are written on the Smooth Tape for input to the 
Profile Program. In Subroutine OZDENS of the Profile 
Program, the standard deviation of In I, sdlni, is computed 
using the formula, sdlni * (sigmaa* + ( sigmab(H-Hb) ) 1 ) * . 
The error associated with the difference of In I at the top 
and the base of the 2-km layer under consideration is then 
the rms sum of the errors associated with each In I. 

The third group of errors are slow drifts that may be 
appreciable during the 15 to 30 minutes required for the 
Rocoz photometer to drop from its ejection altitude to near 
the tropopause (or over the horizon such that telemetry 
signals are no longer received). This Includes changes in 
photodiode responsivity due to changes in temperature or 
length of time in operation, drifts in the zero level of the 
electronics, or slow changes in the absolute calibration of 
the a/d converter. One way of detecting some of these 
changes is to compare the ozone derived from the signal from 
a particular filter to the ozone derived from the ratio of 
that filter signal (counts) to the signal from the SO 
filter, which is at a wavelength for which the ozone 
absorption coefficient is relatively small. In Subroutine 
OVERLP of the Profile Program, a comparison is made between 
the two profiles. 


Altitude 

The error associated with the radar data from the Wallops 
Flight Facility is stated to be 60±30 feet. In Subroutine 
OZDENS of the Profile Program, the error in altitude is set 
to 20 m; this error, expressed in percent, is placed in the 
type 2xx. data record. The error in the vertical thickness 
of the layer under consideration is set to 14 m. 


Air Temperature 

Air temperature is obtained from the Datasonde at altitudes 
above about 20 km and from the ECC sonde at the lower 
altitudes. The error in the Datasonde temperatures is 
typically one degree but increases rapidly above 55 km, 
reaching 15 degrees at 70 km. To this must be added the 
uncertainty due to atmospheric variability. The formulae 
used for this is given in the earlier section on Datasonde 
Observations. The result, expressed in percent, is placed 
in the type 2xx. data record. 
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At altitudes below about 20 km, air temperature is obtained 
from the meteorological sonde to which the ECC sensor is 
attached. The one-sigma error is assumed to be one kelvin. 


Air Pressure 

Air pressure is obtained from the meteorological sonde to 
which the ECC sensor is attached and is assumed to have an 
error of 0.7 mbar. At altitudes above about 20 km, 
pressures are derived from a "base pressure" and the 
datasonde temperature profile, as described in the earlier 
section on Datasonde observations. The error in the base 
pressure is that of the ECC sonde pressure plus errors due 
to atmospheric variability and in matching altitude levels; 
the error in base pressure is assumed to be 1 mbar. At 
higher altitudes the effect of temperature errors must be 
added to this. These estimates are made in Subroutine 
OZDENS, expressed as percent, and placed in the type 2xx. 
data records. 


Ozone Density 

In Subroutine OZDENS the error in the ozone density is set 
to the rms sum of the error in d(ln I) and the error in 
layer thickness, all expressed in percent and the value 
placed in the type 2xx. data record. 

In Subroutine COMPST, the profiles from the individual 
filters are combined to provide a composite profile. At 
each altitude level the available values from the four 
filters are weighted inversely by the associated error. 
That is. 

Composite error « sum( 1 . /err ( i ) ) and 
Composite density * sum (dens i ty (i ) /err (i ) 1 

sum(err ( i ) } 

where i refers to the values available at each level. The 
resulting composite density and error are placed in the type 
300. data records. 


Ozone Overburden 

For the Composite data set, in Subroutine COMPST, the 
initial overburden is obtained by adjusting the column 
content given in the Ozone Model to the density at that 
level in the Composite profile. That is, at the top level, 
Col(comp) « Col(model) * Density (comp) / Density (model) 
The error associated with this is the rms sum of the 
Composite density error and the estimated error in the scale 
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height of the model; this latter error Is assumed to be 
15X. The column content at the next lower level is obtained 
by adding the average density in that layer to the column at 
the top of the layer. The error at this level then is the 
rms sum of the error (in atm-cm) of the column at the top of 
the layer and the errors in the densities in the top half 
and bottom half of the layer. The column error is then 
converted to percent and stored in the type 300. data 
record. 


Ozone Partial Pressure and Mixing Ratio 

In Subroutine CHECKS the Composite data set is completed by 
computing the ozone partial pressure and volume mixing 
ratio. The error associated with the partial pressure is 
the rms sum of the percent errors in the ozone density and 
the air temperature. The error associated with the volume 
mixing ratio is the rms sum of the percent errors associated 
with the ozone partial pressure and the air pressure. The 
Composite data set is then transferred to the type 300. data 
records . 
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Common Areas 


The Profile Program uses seven named common areas to 
facilitate sharing of data among the various subroutines. 
The contents of each common area are described In this 
section. All variables are real (R*4) or Integer (1*4) 
according to the usual naming convention, except where 
specifically Indicated otherwise. The seven named common 
areas are SETUP, PROFLS, FSHAPE , CONSTT, OATAIN, COEFTS, and 
DATOUT. 


COMMON/SETUP/TAPEIN,TAPOUT,NFIN,NFOT,INCLD(7),TOP(7) , 
BASE ( 7 ) 


TAPEIN 

R*8 

Input 

Tape Volume Identifier 

TAPOUT 

R*8 

Output 

Tape Volume Identifier 

NFIN 


File Number on 

input tape 

NFOT 


File Number on 

output tape 

INCLD 

(1) 

03 

process 

S3 filter 


(2) 

02 

process 

S2 filter 


(3) 

01 

process 

SI filter 


(4) 

00 

process 

SO filter 


(5) 

30 

process 

S3/S0 data 


(6) 

20 

process 

S2/S0 data 


(7) 

10 

process 

SI/SO data 


any negative number - ignore this filter 
TOP highest altitude to be processed, km 

BASE lowest altitude to be processed, km 
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COMMON/ PROFLS/DSB (3) ,MXSND,NECC(2) ,ECCLM(3,2) .SONDE ( 5 ,70) 
ECC(4 ,19,2) , ST D ( 3 ,60) 


OSB 


MXSND 

NECC 

ECCLM 


SONDE 


ECC 


STD 


(1) base altitude, km, for use with Datasonde 

(2) temperature (kelvin) at base altitude 

(3) air pressure (mbar) at base altitude 
index of lowest altitude in SONDE 

number of rows in ECC for each of 2 sonr.is 
(1 ,x) ozone column (atm-cm) up to 250 mbar 

(2,x) ozone column (atm-cm) up to 125 mbar 

( 3 ,x) ozone column (atm-cm) up to 62.5 mbar 

for each altitude, starting at highest alt 
(1 ,x) altitude, km, for Datasonde 
(2,x) air temperature, kelvin 
(3,x) air pressure, mbar, hypsometric value , 
(4,x) air density, derived from T and P, kg/nr 
( 5 , x ) error, one sigma, in air temperature, k 

data at standard pressure levels from ECCs 
beginning at 1000 mbar 
(l,x,y) altitude, km 
(2,x,y) air temperature, kelvin 
(3,x,y) air pressure, mbar 
(4,x,y) ozone partial pressure, nanobar 
values from model ozone profile 
(1 ,x) altitude, km, starting at 70 km , 

(2,x) ozone number density, molecules/nr 
(3,x) ozone overburden, atm-cm 


COMMON/ FSHA PE / CW ( 7) ,BW( 7) .CMCBRD ( 4 ,7 ) 


CW 


center wavelength. 

BW 


bandwidth at half- 

CMCBRD 

(1.x) 

AO as given by NWC 


(2,x) 

A 1 


( 3 , x ) 

A2 


(4,x) 

Beta 
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COMMON/CON ST T/AV0GAD,V0LSTP,RSTAR,DGT0RD,GRV45 ,AIRM,GRV, 
RZRO 

AVOGAD 6.022169 26 molecules, Avogadro's number 

per mole 

VOLSTP 22.4136 liters per mole 

RSTAR 8314.32 N m/(kmol K ) , universal gas const. 

DGTORO 1.745329* 2 radians per degree 

GRV45 9.806160 m/s* acceleration of gravity 

at 45 degrees latitude 

AIRM 28.964 kg/kmol, molecular weight for air 

6RV m/s* acceleration of gravity at latitude 

of launch site 

RZRO km, radius of earth at latitude of launch 


COMMON/ DATA I N/ NLA YRS ,SECSZA{60) , SMOOTH ( 20 , 60, 7) 

NLAYRS number of altitude layers used in SMOOTH 

SECSZA secant/Chapman-function of solar zenith 

angle for each layer 

SMOOTH all of the data records from the Smooth 

tape, separated by filter 
Note that in Subroutine OZDENS, the Smooth 
records are replaced by the Profile 
data records for each filter 


COMMON /COE FTS / GAMM { 1000) ,BETA(1000) , ABSZ ( 1000 ) ,ABST(800,2) 

GAMM product of solar flux, diffuser trans- 

mission, and detector responsivity from 
2401 to 3400 angstroms 

BETA Rayleigh scattering coefficient, 2401-3400 

ABSZ ozone absorption coefficient, 2401-3400 A 

at room temperature 

ABST ozone absorption coefficient, 2601-3400 A 

at two lower temperatures 
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COMMON/ DATOUT/NSAV, PE D(20, 30), SAVE(20, 60) 

NSAV The number of records in PED or SAVE to be 

written to the Profile output tape 
PED All the Pedigree records read from the 

Smooth tape or created in Profile 
SAVE Data records created by Profile Program 

In Subroutine COMPST, this area is reassigned thus: 

COMMON/ DATOUT/ NSAV, NDXTOP.NDXBAS.NLVLS , DM ( 1 7 ) .WORK ( 20, 89) 

In Subroutines RELATD and PRTWRT, this common area is 
described as follows: 

COMMON /DATOUT /NSAV, NDXTOP.NDXBAS.NLVLS, DM ( 10 ) ,D0BTM, DO BAM, 
DOBPM , PMDL ( 2 ) .ECTIE ,ECC0L ,W0RK(20 ,89) 


index in WORK or COMP (see below) for top 
of Rocoz ozone profile, 
index for lowest altitude of Rocoz profile 

Rocoz launch date, MMDDYY. 

Rocoz launch time, HHMMSS.S 
launch latitude, ±degrees 
launch longitude, degrees EAST 
Rocoz flight number xxx.x 

launch site, EBCDIC characters 
DM ( 1 1 ) DOBTM date of Dobson observations, MMDDYY. 

DM(12) DOBAM Dobson data, morning, atm-cm 

DM ( 13) DOBPM Dobson data, evening, atm-cm 

DM(14 ,15) PMDL Model identification, EBCDIC characters 

DM ( 1 6 ) ECTIE tie-on altitude for ECC data, km 

DM ( 1 7 ) ECCOL ECC ozone column up to ECTIE, atm-cm 

WORK work space used in preparation of Composite 

ozone profile 

In Subroutines RELATD and PRTWRT, WORK is further 
divided thus 

COMMON/ DATOUT/NSAV ,C0MP(20,70) ,TYP4(20,2) ,TYP5 ( 10,34) 

COMP Rocoz composite ozone profile (type 300. 

data record) 

TYP4 ozone content in std pressure layers 

(type 400. data record) 

TYP5 ozone values at standard pressure levels 

(type 500. data record) 


NSAV 
NDX', OP 

NDXBAS 
NLEVLS 
DM ( 1 ) 

DM ( 2 ) 

DM ( 3) 

DM ( 4 ) 

DM ( 5 ) 

DM ( 6 ) 

DM ( 7 - 1 0 ) 
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Input Data 


The Profile Program uses the Smooth Tape as the source of 
data from the Rocoz photometer; the detailed description of 
this tape is given in the section on Tape Formats. DATA5 
"cards" are the source of specific run options, calibration 
information, and observations from the Dobson 
spectrophotometer, the Datasonde, and ECC sondes. Data 
needed to compute effective absorption and scattering 
coefficients are stored on disc: solar flux, rayleigh 
scattering coefficients, ozone absorption coefficients at 
three temperatures, diffuser transmission, detector 
responsi vity , and, for filters calibrated at the Wallops 
Flight Facility, filter shapes. 


Card Input 


For ease in identification, each card has a letter in the 
first column; this is not read by the Program. 
Approximately 100 cards are needed, depending upon the 
altitude range of the Datasonde and the number of ECC 
sondes . 


Columns Format Content Variable 


1 A One card required. Tape and Filter IDs 


2-7 

F6.1 

Flight number and segment 

FLT(l) 

10-17 

A8 

DT, Tape Volume Identifier, 
volume serial number for 
input (Smooth) tape 

TAPEIN 

18-20 

13 

file number on input tape 

NFIN 

25-32 

A8 

DT, Tape Volume Identifier, 
output (Profile) tape 

TAP0UT 

33-35 

13 

file number on output tape 

NF0T 

40-41 

12 

month, filter calibration 

MDY(l) 

43-44 

12 

day, filter calibration 

MDY( 2 ) 

46-47 

12 

day of month, filter calibration 

MDY ( 3 ) 

49-52 

A4 

filter manufacturer 

MFG 

54-58 

15 

dddyy, date of filter set 
assembly 

ICAL 

60-65 

16 

filter set identification no. 

ICLID 
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Columns 

1 

3-4 

5-9 

10-14 

17-22 

24-28 

31-38 

40-49 

51-60 

62-68 

Notes : 


1 

3-4 

6-7 

9-10 

13-17 

24-24 


Card Input (continued) 


Format Content Variable 


B Seven required. Filter processing 


12 

03, 02, 01, 00, 30, 20, or 10 

; INCLD(x) 


use negative number (-1) to 

skip 

F5.1 

highest altitude to be used 

T0P(x) 

F5.1 

lowest altitude to be used 

BASE ( x ) 

F6.2 

center wavelength, nm 

CW( x ) 

F5.2 

bandwidth, nm 

BW ( x ) 

F8.4 

McBride's A0 

CMCBRD ( 1 ,x) 

F10.5 

McBride's Al 

CMCBRD ( 2 , x ) 

F10.5 

McBride ' s A2 

CMCBRD ( 3 , x ) 

F7.5 

McBride's B 

CMCBRD ( 4 , x ) 


- Values in INCLD must be in the order indicated. 

- For cards 5, 6, and 7, CW and BW are blanks. 

- The subscript, x, refers to card number (1-7) 

- The values of center wavelength, bandwidth, AO, 
Al, and A2 are also calculated for the filters 
calibrated at WFF; B is obtained for the center 
wavelength from the Rayleigh scattering table. 


C One required. Dobson data 

12 month, date of Dobson obs. MDY(l) 
12 day of month, date of Dobson obs. MD Y ( 2 ) 
12.1 year , date of Dobson observations M D Y ( 3 ) 
F5.1 AM value, in Dobson units PE D ( 3 , 6 ) 
F5.1 PM value, in Dobson units PED(4,6) 
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Columns 

1 

3-10 

13-14 

16-17 

19-20 

23-24 

25-26 

29-31 

34-40 

43-48 

51-56 

Note: 


1 

3-8 

11-16 

19-22 

Notes : 


Card Input (continued) 


Format Content Variable 


0 One required. Datasonde ID 

A8 Datasonde flight number. FLTID 

12 month, date of datasonde flight MDY(l) 

12 day of month, date of datasonde MDY ( 2 ) 

12 year, date of datasonde MDY ( 3 ) 

12 hour, time of datasonde launch IH 

12 minute, time of datasonde launch IM 

13 minutes between datasonde and Rocoz MDEL 

launches 

F7.1 base altitude, meters DSB(l) 

F6.2 base temperature, C. or kelvin DS B ( 2 ) 

F6.2 base pressure, mbar DSB(3) 


Base altitude, temperature, and pressure are obtained 
from the study of data from the nearly simultaneous 
ECC sondes and are used to calculate the pressures 
and densities at Datasonde altitudes. 


E Max of 70 required. Datasonde Data 

F6.0 altitude of observation, m or km ALT 

F6.2 observed temperature, C. or kelvin TEMP 

F4.1 error in temperature, 2-sigma, deg TERR 

- First card gives data for the highest altitude. 

- Last card must set "altitude" to -1, to indicate 
end of data. 

- If the Datasonde report does not include error 
estimate, use errors from F.Schmidlin (JGR.86, 
9599,1981) 


11-71 


PROFILE PROGRAM 


Card Input (continued) 


Columns Format Content Variable 


1 


F One 

required. ECC sonde ID 


3-10 

A8 

ECC sonde flight number. 

FLTID 

13-14 

12 

month , 

date of ECC sonde flight 

MDY(l) 

16-17 

12 

day cf 

month, date of ECC sonde 

MDY ( 2 ) 

19-20 

12 

year. 

date of ECC sonde 

MDY (3) 

23-24 

12 

hour , 

time of ECC sonde launch 

IH 

25-26 

12 

minute 

, time of ECC sonde launch 

IM 

29-34 

F6.0 

tie-on 

altitude, meters 

TIEALT 

37-43 

F7.5 

ozone 

column up to TIEALT, atm-cm 

TIEC0L 


Note: If data from a second ECC sonde is available, a type 

F card for the second sonde follows the type G 
(data) cards for the first sonde. 


1 G 19 required. ECC sonde Data 


3-7 

15 

altitude, meters 

KH 

10-14 

F5.1 

ozone pressure, nanobars 

ECC ( 4 , 1 , 1 ) 

17-22 

F6.1 

air pressure, millibars 

ECC ( 3 , I , 1 ) 

25-29 

F5.1 

air temperature, C. or kelvin 

ECC ( 2 , I , 1 ) 


Note: Data are given at 19 standard pressure levels, namely, 

1000., 850., 700., 500., 400., 300., 250., 200., 150., 

100., 70., 60., 50., 40., 30., 25., 20., 15., and 10. 
If data are not available for a given pressure level, 
the altitude, ozone pressure, and temperature are to 
be set to 0.0 The first card is for 1000. mbar. 


If data from a second ECC sonde are available, they 
follow the type F card identifying the flight. 


1 


H One required. ECC Ozone Column 


7-13 

20-26 

33-39 


F7.5 

ECC 

03 

column , 

up 

to 

250 

mbar 

ECCLM(l) 

F7.5 

ECC 

03 

column , 

up 

to 

125 

mbar 

ECCLM ( 2 ) 

F7.5 

ECC 

03 

column , 

up 

to 

62. 

5 mbar 

ECCLM ( 3) 


Note: Column content is in atm-cm of ozone. 
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RADAR TAPE FORMATS 1 1 1 -2 

Radar Tape III-2 

Radar (1983) Tape 1 1 1 -4 

REFORMATTED RADAR TAPE 1 1 1 -5 

TELEMETRY TAPE III-7 

EDITED OZONE DATA TAPE 1 1 1 -10 

MERGE TAPE 1 1 1 - 1 2 

SMOOTH TAPE 1 1 1 -14 

PROFILE TAPE I X I -16 

Pedigree Records 1 1 1 - 1 8 

Data Records 1 11-23 
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RADAR TAPE FORMAT 


Most tapes for flights through April 1982 (Rocoz flight no. 
315) were prepared by the Wallops Station Meteorological 
Support program 1.1.2123. The "GE-625 Meteorological 
Support Format" (MESUP) is as follows: 

Each physical record is composed of twenty-one 
twel ve-character words. The word layout is in the 
following form: 

XXXXXXXXXX S 0 
0123456789 10 11 


Character number 0 is the most significant character and 
character number 11 is the least significant character 
and will always be a zero. The sign position is fixed 
for all words and is in character position 10. The 
first word of each record is a system word recorded for 
tape identification and may be ignored in all 
processing. 

The tape is recorded with word one character number 0 
being recorded first. Each tape file will be followed 
by a tape mark. Inter-record gap is approximately 
three-fourths of an inch. Each tape will have a ten to 
twenty foot leader of blank tape before the beginning of 
tape (BOT) tinsel marker. 


Tape Description and Parameters for JCL 


7-track 
no label 

records are fixed length 
block size is 252 bytes 
record length is 252 bytes 
density * 800 bpi 
even parity 
BCD characters 


UNIT*2400-7 or 7TRACK 

LABEL=NL 

RECFM=F 

BLKSI ZE =252 

LRECL=252 

DEN = 2 

TRTCH ®E 

TRTCH-T (to EBCDIC) 
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Radar Tape Format 


Word 

Character 

Decimal point 

assumed 

Number 

Positions 

Data Content between characters 

2 

0-9 

Elapsed time given in seconds. 

5&6 

3 

0-9 

Slant range from the origin in feet. 

9 & 1 0 

4 

0-9 

Azimuth from origin given in degrees. 

5 & 6 

5 

0-9 

Elevation from origin given in degrees. 

5 & 6 

6 

0-9 

Horizontal range from origin in feet. 

9&10 

7 

0-9 

North-South range from origin in feet. 

9&10 

8 

0-9 

East-West range from oriT*n in feet 

9&10 

9 


(Not used, contains zeros) 


10 

0-9 

Altitude above an oblate spheroid in feet. 

9&1 0 

11 


(winds or zeros) 


12 


(winds or zeros) 


13 

0-9 

Horizontal component of earth relative 

9&10 



velocity vector in feet/second .* 


14 

0-9 

East-West component of negative velocity 

9&10 



vector in feet/second.* 


15 

0-9 

North-South component of negative velocity 

9&10 



vector in feet/second.* 


16 

0-9 

Altitude component of velocity vector 

9&10 



in feet/second.* 


17 


(Not used, contains zeros) 


18 

0-9 

Azimuth of the horizontal component of the 

5&6 



negative velocity vector in degrees.* 


19 


(Not used, contains zeros) 


20 


(Not used, contains zeros) 


21 


(not used, contains zeros) 



* Velocity vector is defined as positive in the direction of 
target motion measured clockwise from North. 
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RADAR (1983) TAPE 


Radar data for flights in 1983 and later are prepared on the 
Applications Directorate Eclipse computer at WFF and provide 
location in meters, latitude and longitude on a 9-track 
tape. This PASS1 tape may be listed by the Listrad Program. 
The Merge Program has been modified such that it can accept 
either this tape or the Reformatted Radar Tape produced by 
the Radar Program. When this, the Radar (1983) tape, is 
used, the Merge Program prepares the Radar Pedigree (Type 
-100.) record. 

Tape Description and Parameters for JCL 


9-track 
no label 

records are fixed length 
and blocked 
blocksize is 32 bytes 
record length is 32 bytes 
density = 800 bp i 


UNIT-800 
LABE L-NL 
RECFH-FB 

BLKSI ZE-32 

LRECL*32 

DEN-2 


Radar (1983) Tape Format 
All words are binary double precision, REAL*8 . 
Word Number Data Content 

1 Seconds since launch time 

2 Altitude, meters 

3 Latitude, degrees North 

4 Longitude, degrees East 
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REFORMATTED RADAR TAPE 


The Reformatted Radar Tape is a product of the Radar 
program. Like all tapes produced by this Rocoz set of 
programs, the first record is a "Pedigree" record giving the 
flight number, date of creation of this tape, and other 
Rocoz-related data entered via the DATA5 input unit. The 
program converts times to seconds after the hour instead of 
after launch and converts English units to metric units. 
Each file contains all available data for a single flight. 


Tape Description and Parameters for JCL 


9 track 

standard label 
records are fixed length 
and blocked 

blocksize is 5200 bytes 
record length is 52 bytes 
tape density is 1600 bp i 


UNIT-2400-9 or 9TRACK 
LABEL-SL 

RECFM-FB 
BLKSI ZE-5200 
LRECL-52 
DEN = 3 


Pedigree Record 


Words 1-9 are R*4 real words; 10-13 are EBCDIC characters. 


Word Content 


Source 


1 -100. Assigned 

2 Flight Number DATA5 

3 VXX . where V indicates vehicle type: DATA5 

*1 rocket 
-2 balloon 

and XX indicates launch site: 

-01 Wallops Island, Virginia 
-02 Poker Flats, Alaska 
*03 Natal , Brazil 
*04 Marambio, Antarctica 
-05 Palestine, Texas 


4 MMDDYY . month, day, year of launch DATA5 

5 HHMMSS.S hour, minute, second of launch DATA5 

6 ±XX.XXX latitude of launch site DATA5 

7 XXX. XXX longitude of launch, degrees WEST DATA5 

8 MMDDYY. date radar tape was received DATA5 

9 MMDDYY. date of creation of this tape Computer 

10-13 launch site, 16 EBCDIC characters DATA5 
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Data Records (Reformatted Radar Tape) 
All words are R*4 real words 


Word Content 


Unit 


1 time from hour prior to start of data seconds 

2 slant range meters 

3 azimuth from origin degrees 

4 elevation from origin degrees 

5 horizontal range from origin meters 

6 north-south range from origin meters 

7 east-west range from origin meters 

8 altitude above oblate spheroid meters 

9 horizontal component of velocity relative 

to the earth meters /second 

10 east-west component of velocity meters /second 

11 north-south component of velocity meters /second 

12 vertical velocity meters/second 

13 azimuth of horizontal component of velocity degrees 


The velocity vector is defined as positive in the direction 
of target motion measured clockwise from North. 
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TELEMETRY TAPE FORMAT 


The telemetry tapes for the Rocoz flights have been produced 
by the "Digital Telemetry System (PDP 11/60): DITES" at the 
Wallops Flight Facility. NOTE THAT THE BYTES FROM THE 
PDP-11 MACHINE ARE PACKED MSB LAST. 


Tape Description and Parameters for JCL 


9-track 
no label 

record length is undefined 
block size is 2440 bytes 
density is 800 bpi 


UNIT *2400-9 or 9TRACK 

LABEL-NL 

RECFM-U 

BLKSI ZE *2440 

DEN-2 


The "televent-ll" data tapes from this WFF system consist of 
three different types of records. The first record of the 
Rocoz data tape is a Header Record that contains up to 80 
ASCII characters. The second type of record is a Data 
Description Record which describes the data for the Rocoz 
data scream. The third type of record is the Data Record . 


Header Record 


Word 

Number Content Description 


1 177776 

2 NN 

3 000000 

4 B A 

5 D C 

6 N 

7 B A 

through 

46 2 Y 


-2 for DEC Fortran compatibility 
number of words In record 
0 to identify Header Record 
4-character stream Identification 
(note order or characters) 
tape number 

ASCII text. Each character in a 
word is recorded most significant 
bit last. 
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Data Description Record 


Word Typical 

No . Content Val ue Descri ption 


1 

177776 

-2 

-2 for DEC Fortran compatibility 

2 

N 

20 

number of words in record 

3 

177775 

-3 

-3 to identify Data Description Record 

4 

B A 

PC 

4-character Run Identifier from setup 

5 

D C 

M4 

of tape processing 

6 

B A 

PC 

4-character Stream Identifier (usually 

7 

D C 

M4 

the same as words 4 and 5 

8 

N 

8 

number of data words per frame 

9 

N 

no 

number of frames per buffer 

10 

N 

0 

scaling factor (normally = 0) 

11 

N 

3 

number of frame Merge words 

12 

N 

11 

number of words/frame 

13 

N 

2 

number of preface words above 8 

14 

N 

0 

buffer appendix size 

15 

N 

1210 

INPUT buffer size (words) 

16 

• i 

1220 

OUTPUT buffer size (words) 

17 

N 

0 

fr ames/subframe number 1 

18 

N 

0 

frames/subframe number 2 

19 

N 

0 

recording mode (0=wc>d; l=byte) 

20 

020040 


two ASCII blanks 


Data Record 


The data record consists of a 10-word preface, 110 11-word 
data frames, and no appendix. 

Preface 


Word Typical 

No . Content Value Description 


1 

17 7776 


2 


N 

12:* 

3 


N 


4 

B 

A 

PC 

5 

D 

C 

M4 

6 

A 

P 

10 

7 

L 

M 

1210 

O 

S 

F 


9 


N 

no 

10 


N 

0 


-2 for DEC Fortran compatibility 
record length (words) 
block count 

4-character Stream Identifier 

A is no. of Appendix words ( ); 

P is no of Preface words ( ) 
INPUT buffer size 
TFE status and frame information 
(see below for explanation) 
number of frames per buffer 
scaling factor (normally * 0) 
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Word no. 8: TFE status 


B i t Content 

0-8 frame number of first search 
9 sub-corn 2 search on last buffer 

10 sub-corn 1 search on last buffer 

11 frame search on last buffer 

12 multiple search on one of below 

13 sub-corn 2 search flag 

14 sub-corn 1 search flag 

15 frame search flag (most significant bit) 

A hex value of 00 is normal; FF usually indicates missing 
frames . 


Word Number 


1 
2 

3 

4 

5 

6 

7 

8 
9 

10 
11 

Word 9: 

: 1 5 14 13 12 111 10 

'tenths of sec. {hundredths of s {thousandths s { 

Word 10: 

15 {14 13 12 111 10 9 8 { 7 6 5 4 { 3 2 1 0{ 

{tens of min{ units of min { tens of sec {units of sec.' 

Word 11: 

{15 14 {13 12 11 10 { 9 8 7 6 { 5 4 { 3 2 1 0{ 

JlQO's {tens of day { units of day {tens hrjunits of hour} 


Data Frame 

Content 

S3 filter, counts 

S2 filter, counts 

compensation filter, counts 

SI filter, counts 

SO filter, counts 

battery voltage, counts 

marker pulse, counts 

uncompensated value, counts 

milliseconds (see telow) 

minutes and seconds (see below) 

day of year, hour of day ( GMT ) (see below) 


9 8 {7 6 5 4 { 3 2 10 
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EDITED OZONE DATA TAPE FORMAT 


The Edited Ozone Data Tape is the product of the Edit 
Program. Following the single Pedigree record are Data 
records, one record for each rotation of the filter wheel. 
Battery voltage and temperature data are converted to 
engineering units; the remaining data are in telemetry 
counts. All words are R*4 real numbers. 


Tape Description and Parameters for JCL 


9 track 

standard label 
records are fixed length 
and blocked 

blocksize is 8000 bytes 
record length is 80 bytes 
tape density is 1600 bpi 


UNIT=2400-9 or 9TRACK 
LABEL=SL 

RE CFM = FB 
BLKS I ZE =8000 
LRECL=80 
DE N = 3 


Edit Pedigree Record 


Word 

Content 



Source 

1 

-200 . 



Assigned 

2 

XXX. s 

flight number 

. segment 

DATA5 



for rocket flights, S is 0 


3 

Payload 

number 


DATA5 

4 

MMDDYY . 

month, day. 

year that telemetry 

DATA5 


tape 

was received. 



5 

MMDDYY. 

month, day. 

year for creation of 

Computer 


th i s 

tape . 



6 

Battery 

cal ibration : 

volts (upper level) 

DATA5 

7 

M 

ll 

tm count *' " 

DATA5 

8 

H 

II 

volts (lower level) 

DATA5 

9 

II 

ll 

tm count " " 

DATA5 

10 

Temperature calib.: 

degree C. (lower) 

DATA5 

11 

N 

•* 

tm count " 

DATA5 

12 

II 

•1 

degree C. (mid) 

DATA5 

13 

N 

It 

tm count " 

DATA5 

14 

II 

II 

degree C. (high) 

DATA5 

15 

H 

n 

tm count 

DATA5 

16 

0.0 




17 

0.0 




18 

0.0 




19 

0.0 




20 

0.0 
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Edit Data Record 


Word 

Content 



Unit 

1 

time since hour of tm word 

1 (SO) 


seconds 

2 

time since hour of tm word 

2 (SI) 


seconds 

3 

time since hour of tm word 

4 ( S2 ) 


seconds 

4 

time since hour of tm word 

5 (S3) 


seconds 

5 

value of tm word 1 (SO) 



counts 

6 

value of tm word 2 (SI) 



counts 

7 

value of tm word 4 (S2) 



counts 

8 

value of tm word 5 (S3) 



counts 

9 

value of marker pulse (tm i 

word 7) 


counts 

10 

value of base level (tm word 71 


counts 

11 

compensation word (tm word 

3) for 

SO 

counts 

12 

compensation word (tm word 

3) for 

SI 

counts 

13 

compensation word (tm word 

3) for 

S2 

counts 

14 

compensation word (tm word 

3) for 

S3 

counts 

15 

length of this cycle 



seconds 

16 

average length of last ten 

cycles 


seconds 

17 

average value of word 3 



counts 

18 

battery voltage (average over 110 

frames ) 

volts 

19 

temperature (from marker pulse height) 

degrees C. 

20 

number of frames in this cycle 


count 
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TAPE FORMATS 


MERGE TAPE FORMAT 


The Merge Tape is the product of the Merge Program. The 
first record is the Pedigree record from the Reformatted 
Radar Tape; the second record is the Pedigree record from 
the Edited Ozone Tape. The third record is the Pedigree 
record created by the running of the Merge Program. The 
data records essentially replace the time information with 
altitude information. All words on this tape are R*4 real 
numbers . 


Tape Description and Parameters for JCL 


9 

track 

UNIT=2400-9 or 

9TRACK 

standard label 

LABE L = SL 


records are fixed length 




and blocked 

RECFM=FB 


blocksize is 8000 bytes 

BLKSI ZE *8000 


record length is 80 bytes 

LRE CL=80 


tape density is 1600 bpi 

DEN = 3 



Merge Pedigree Record 


Word 

Content 


Source 

1 

-300 . 


Assigned 

2 

XXX. S flight number . segment 

DATA5 


for rocket flights, S is 0 


3 

MMDDYY. month, day, year 

for creation of 

Computer 


this tape. 



4 

Apparent right ascension 

of the sun, hour 

DATA5 

5 

ii II ll 

minutes 

DATA5 

6 

7 

If ll ll 

ii ll ll 

seconds 

DATA5 

/ 

difference between two 

t 

days, seconds 

DATA5 

8 

Apparent declination of 

the sun , degrees 

DATA5 

9 

li ll 

, minutes 

DATA5 

10 

U N 

, seconds 

DATA5 

11 

N II 

• 



difference between two 

days, seconds 

DATA5 

12 

Apparent sidereal time. 

hour 

DATA5 

13 

«i ll 

minutes 

0ATA5 

14 

N II 

seconds 

DATA5 

15 

0.0 



16 

0.0 



17 

0.0 



18 

0.0 



19 

0.0 



20 

0.0 




1 1 1 -12 



TAPE FORMATS 


Merge Data Record 







Word 






on Edit 

Word 

Content 



Unit 

Data Rec 

1 

time since hour of tm word 1 

(SO) 


seconds 

i 

2 

length of this cycle 



seconds 

15 

3 

number of samples in this cycle 


count 

20 

4 

temperature of instrument 



degrees 

C. 19 

5 

value of tm word 1 (SO) 



counts 

5 

6 

value of tm word 2 (SI) 



counts 

6 

7 

value of tm word 4 (S2) 



counts 

7 

8 

value of tm word 5 (S3) 



counts 

8 

9 

value of marker pulse (tm word 7) 


counts 

9 

10 

battery voltage 



volts 

18 

11 

compensation word (tm word 3) 

for 

SO 

counts 

11 

12 

compensation word (tm word 3) 

for 

SI 

counts 

12 

13 

compensation word (tm word 3) 

for 

S2 

counts 

13 

14 

compensation word (tm word 3) 

for 

S3 

counts 

14 

15 

average value of compensation 

word 


counts 

17 

16 

altitude for tm word 1 (SO) 



meters 


17 

altitude for tm word 2 (SI) 



meters 


18 

altitude for tm word 4 ( S 2 ) 



meters 


19 

altitude for tm word 5 (S3) 



meters 


20 

solar zenith angle 



degrees 
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SMOOTH TAPE FORMAT 


The Smooth Tape is the product of the Smooth Program. The 
first record is the Pedigree record from the Reformatted 
Radar Tape; the second record is the Pedigree record from 
the Edited Ozone Tape. The third record is the Pedigree 
record from the Merge Program. The fourth record is the 
Pedigree record created by the Smooth program. Each Smooth 
Data record provides data at one integer altitude for one 
filter. After data records for each of the four filter?, 
similar information is given for S3/S0, S2/S0, and SI/SO. 
All words are R*4 real numbers. 


Tape Description and Parameters for JCL 


9 track 

standard label 
records are fixed length 
and blocked 

blocksize is 4000 bytes 
record length is 80 bytes 
tape density is 1600 bpi 


UNIT=2400-9 or 9TRACK 
LABEL=SL 

RECFM=FB 
BLKSI ZE*4000 
LRECL-80 
DEN=3 


Smooth Pedigree Record 


Word Content 


Source 


1 -400. Assigned 

2 XXX. S flight number . segment DATA5 

for rocket flights, S is 0 

3 MMDDYY. month, day, year for creation of Computer 


this tape. 


4 

zero offset for SO, counts 



DATA5 

5 

zero offset for SI , counts 



DATA5 

6 

zero offset for S2 at temp. 

of 

word 11 

DATA5 

7 

zero offset for S3 at temp. 

of 

word 11 

DATA5 

8 

maximum integer altitude to 

be 

processed 

DATA5 

9 

minimum integer altitude to 

be 

processed 

DATA5 

10 

minimum acceptable value of 

comp, word 

DATA5 

11 

temperature (deg. C) for words 

6 and 7 

DATA5 

12 

temperature (deg. C) for words 

13 and 14 

DATA5 

13 

zero offstt for S2 at temp. 

of 

word 12 

DATA5 

14 

zero offset for S3 at temp. 

of 

word 12 

DATA5 

15 

0.0 




16 

0.0 




17 

0.0 




18 

0.0 




19 

0.0 




20 

0.0 
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Smooth Data Record 


Word Content 


Unit 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


filter position: 3., 2., 1., or 0. 

if ratio to SO: 30., 20., or 10. 
altitude level, H 

signal at XCEN, I=EXP(A+B*(H-Hb) ) . 

error = (SIGMAA 1 * ( H-Hb) * S I GMAB * ) * 
solar zenith angle at XCEN 
time at XCEN 

temperature of instrument 

average length of filter cycle (99 cycles) 

battery voltage 

average compensation word (99 cycles) 
upper end of altitude interval 
lower end of altitude interval 
number of points used 

estimated signal at Hb altitude, *EXP(A) 
standard deviation of A (SIGMAA) 
estimated absorption and scatter i ng ,=-B 
standard deviation of B (SIGMAB) 
standard deviation of fit to curve 
correlation coefficient 

number of points provided in altitude interval 


ki 1 ometers 
count or ratio 
See note, 
degrees 

seconds after hour 
degrees C 
seconds 
volts 
counts 
ki 1 ometers 
ki 1 ometers 
count 


Note: In files created before Dec. 1983, word 4 is slope {= B*I). 
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TAPE FORMATS 


PROFILE TAPE FORMAT 


The Profile Tape is the product of the Profile Program and is 
Intended to be suitable for the archival of data from the Rocoz 
flight. It includes not only the derived ozone profile, but 
also the related ECC and Datasonde data and the calibration 
constants used in the processing of the Rocoz data. The first 
group of records are Pedigree records, both those from 
preceding programs in the Rocoz data reduction system and 
created by the Profile program. The second group of records 
give ozone content as a function of altitude with values of 
various parameters developed in the computation. The third 
group of records give the ozone content in various units and 
altitude levels for the convenience of the user. The maximum 
number of records for a flight is 456. Most of the data are in 
R*4 real words except for some words in the Pedigree records 
which are in character (EBCDIC) format. 


Tape Description and Parameters for JCL 


9 track 

standard label 
records are fixed length 
and blocked 

blocksize is 4000 bytes 
record length is 80 bytes 
tape density is 1600 bpi 


UNIT-2400-9 or 9TRACK 
LABEL-SL 

RECFM-FB 
BLKSI ZE-4000 
LRECL-80 
DEN-3 
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Types of Records 


Number of 


First Word 

Content 

this 


-100. 

Pedigree: reformatted radar tape 


1 

-200 . 

Pedigree: edited telemetry tape 


1 

-300. 

Pedigree: merge data tape 


1 

-400. 

Pedigree: smooth data tape 


1 

-500. 

filter: McBride coefficients 


1 

-501. 

filter coeff; Dobson data 


1 

-502 . 

Datasonde flight 10 and data 


1 

-503. 

Datasonde data 


3 

-504 . 

ozone density model 


1 

-505 . 

ozone model, continued 


3 

-506. 

ECC sonde ID and data 


1 

-507 . 

ECC data, continued 


4 

-506. 

second ECC ID and data, if available 


1 

-507 . 

second ECC data, continued 


4 

-508. 

source of data for solar flux 


1 

-509. 

source of data for diffuser transmission 

1 

-510. 

source of data for detector responsivity 

1 

-511. 

source of Raleigh scattering coefficient 

1 

-512. 

source of ozone absorption coeff., max 

temp 

. 1 

-513. 

source of ozone absorption coeff.. 


1 

-514. 

source of ozone absorption coeff., min 

temp 

. 1 


For each filter, at each altitude level, are two aata records; 
the number of records per filter or filter ratio varies from 0 
to 25. 

lxx. intermediate values in the calculation of ozone 
2xx. observed ozone with air temp., air press., etc. 


The following records give derived profiles 
300. composite ozone profile, 

1 km pe>'' record 

400. ozone layer content in 12 std. 
pressure layers 

500. ozone at standard (33) pressure 

2 levels per record. 


in various units: 
10-60 
2 

levels , 

17 


Maximum number of records: 460. 
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Pedigree Records Created by Profile Program 


Word Content 

1 -500. Rocoz flight ID; filter ID end calibration 

2 XXX. S flight number and segment 

3 MMDDYY . month, day, year for creation of this tape 

4 number of ECC sondes available: 0., 1., or 2. 

5 AAAA (EBCDIC) manufacturer of filter set 

6 DDDYY . day of year and year of set assembly 

7 filter set identification number 

8 MMDDYY. month, day, and year of filter calibration 

9 S3 center wavelength, nm 

10 S3 bandwidth, nm 

11 S3 AO, ozone absorption coeff. from McBride's calibration 

12 S3 Al, coeff. for rate of change due to ozone overburden 

13 S3 A2, 

14 S3 beta, Rayleigh scattering coefficient 

15 S2 center wavelength, nm 

16 S2 bandwidth, nm 

17 S2 AO, ozone absorption coeff. from McBride's calibration 

18 S2 Al, coeff. for rate of change due to ozone overburden 

19 S2 A2 , " 

20 S2 beta, Rayleigh scattering coefficient 


1 -501. Dobson data; filter calibration data 

2 MMDDYY. month, day, year of Dobson observations 

3 AM data, Dobson units; 0.0 If unavailable 

4 PM data, Dobson units; 0.0 if unavailable 

5 0.0 

6 0.0 

7 0.0 

8 0.0 

9 SI center wavelength, nm 

10 SI bandwidth, nm 

11 SI A0, ozone absorption coeff. from McBride's calibration 

12 SI Al, coeff. for rate of change due to ozone overburden 

13 SI A2, 

14 SI beta, Rayleigh scattering coefficient 

15 SO center wavelength, nm 

16 SO bandwidth, nm 

17 SO A0, ozone absorption coeff. from McBride's calibration 

18 SO Al , coeff. for rate of change due to ozone overburden 

19 SO A2, 

20 SO beta, Rayleigh scattering coefficient 
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Word Content 


1 

2 

3 

4 

5 

6 
7 
6 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


-502 . Oatasonde ID and data 

initial altitude (top) of datasonde profile, km 

end altitude (bottom) of datasonde profile, km 

number of levels in datasonde profile 

AAAA (EBCDIC) Datasonde flight identification 

AAAA (EBCDIC) N (continued) 

MMDDYY . month, day, year of datasonde flight 

HHMM . hour, minute of datasonde launch 

base altitude for pressure profile, kilometers 

air pressure at base altitude, millibar 

acceleration of gravity at latitude of launch site, m/s* 

earth radius at launch latitude, km 

temperature at initial altitude, kelvin 

temperature at next lower altitude, kelvin 

M 


ll 


1 -503. Datasonde temperature data, continued 

2 temperature at next lower altitude, kelvin 
3-20 

There may be more than one type -503. record, depending on the 
number of levels at which data are available. Where no more 

levels are available, the remainder of the words in the record 
are set to 0.0 
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Word Content 

1 -504. Model values for ozone overburden 

2 initial altitude of model, km (typically, 70.) 

3 altitude increment, km (typically, -1.) 

4 number of levels 

5 AAAA (EBCDIC) model Identification 

6 AAAA (EBCDIC) (continued) 

7 altitude for Initial overburden of ozone, km 

8 ozone overburden above altitude In word 7, number/m* ^ 

9 ozone number density at the initial altitude, number/m 

10 ozone number density at next lower altitude 

11-20 


1 -505. Model values, continued 

2-20 ozone number density at next lower altitude, or padded 
with 0.0 

There may be 0 - 3 type -505. records, depending upon the model 
in use. 


1 -506. ECC flight ID and data 

2 tie-on altitude, typically 20 to 25 km 

3 ozone column content, ground to tie-on altitude, atm-cm 

4 number of levels of ECC data 

5 AAAA (EBCDIC) ECC sonde flight Identification 

6 AAAA (EBCDIC) (continued) 

7 MMDDYY. month, day, year of ECC flight 

8 HHMM. hour, minute of ECC launch 

9 altitude, beginning near the ground, km 

10 air temperature, kelvin 

11 air pressure, mbcr 

12 ozone partial pressure, nanobar 

13-20 repeat parameters of words 9-12 for increasing altitude 
levels . 

Data are to be given at 19 standard pressure levels: 

1000. 850. 700. 500. 400. 300. 250. 200. 150. 

100. 70. 60 50. 40. 30. 25. 20. 15. 10. 

If data from a second ECC sonde are available, record types -506. 
and -507. are to be repeated for the second flight. 
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TAPE FORMATS 



f 


Word Content 

1 -507. ECC data, continued 

2 1. or 2., first or second ECC sonde 

3 AAAA (EBCDIC) ECC sonde flight identification 

4 AAAA (EBCDIC) (continued) 

5-20 repeat (altitude, temperature, pressure, ozone) 

as needed for the remaining standard pressures. 

There may be up to four cards of this type for each ECC sonde. 


1 -508. source of data for solar flux 

2-17 alphanumeric description of data 

18 initial wavelength of data, A. 

19 end wavelength of data, A. 

20 0.0 


1 -509. source of data for diffuser transmission. 

2-17 alphanumeric description of data. 

18 initial wavelength of data, A. 

19 end wavelength of data, A. 

20 0.0 


1 -510. source of data for detector respons i v i ty . 

2-17 alphanumeric description of data 

18 initial wavelength of data, A. 

19 end wavelength of data, A. 

20 0.0 


1 -511. source of data for Rayleigh scattering coefficient 

2-17 alphanumeric description of data 

18 initial wavelength of data, A. 

19 end wavelength of data, A. 

20 0.0 


1 -512. source of data for ozone absorption coefficient, 

room temperature. 

2-17 alphanumeric description of data 

18 initial wavelength of data, A. 

19 end wavelength of data, A. 

20 temperature, K 
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1 '513. source of data for ozone absorption coefficient, 

lower than than in Pedigree record *512. 

2*17 alphanumeric description of data 

18 initial wavelength of data, A. 

19 end wavelength of data, A. 

20 temperature, K 


1 -514. source of data for ozone absorption coefficient, 

lowest available temperature. 

2*17 alphanumeric description of data 

18 initial wavelength of data, A. 

19 end wavelength of data, A. 

20 temperature, K 
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Data Records Created by Profile Program 


Word Content 

1 lOf. for word 5 In counts; lfO. word 5 as ratio to SO 

f * filter number: 3, 2, 1, or 0 

2 h, altitude, km 

3 compensation channel, counts 

4 secant /Chapman -function for solar zenith angle 

5 (I), intensity at altitude, h, counts or ratio to SO 

6 slant ozone overburden, atm-cm 

7 slant air mass overburden, atm 

8 effective ozone absorption coefficient at h, alpha-eff 

9 effective Rayleigh scattering coefficient at h, beta-eff 

10 d ( 1 n I ) / d h 

11 d alpha-eff/d slant-ozone 

12 d slant-air/dh 

13 optical depth due to slant ozone 

14 optical depth due to slant air 

15 Fraser correction (for scattered light) 

16 noise for delta (In I) 

17 delta-h used by Smooth in deriving I(top) - I(base) 

18 B, (absorption ♦ scattering) from Smooth 

19 std. dev. of B 

20 photometer temperature, degree C 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


20f. from this filter only; 2f0. from ratio to SO 
f* filter number: 3, 2, 1, or 0 
h, altitude, km 
error in h, percent 

HHNNSS.S, hours, minutes, seconds, time, GMT 
solar zenith angle, degrees 

noise in ozone density, percent, (from d ( 1 n I ) / dh noise) 

ozone density, atm-cm/km 3 

ozone density, number/nr 

ozone vertical overburden, atm-cm 

ozone vertical overburden, number/m* 

air temperature, kelvin 

error in air temperature, percent 

air pressure, millibar 

error in air pressure, percent 

air density, number/m 3 
air, vertical column, atm. 
ozone, partial pressure, nanobar 
ozone volume mixing ratio 
seconds after the hour 

record number, except for last for this filter, then *-l. 


Note: All errors are for one standard deviation. 
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The following records appear after all lrf. and 2rf. records for 
all filters and ratios have been written. 

Word Content 

1 300. Composite Ozone Table (based on all Rocoz filters) 

2 h , altitude, km 

3 ozone density, atm-cm/km ^ 

4 ozone density, number/m' 1 

5 ozone density, error, percent 

6 ozone vertical column, atm-cm 

7 ozone vertical column, error, percent 

8 air temperature, kelvin 

9 air temperature, error, percent 

10 air pressure, millibar 

11 air pressure, error, percent 

12 ozone partial pressure, nanobar 

13 ozone partial pressure, error, percent 

14 ozone volume mixing ratio, ppm 

15 ozone volume mixing ratio, error, percent 

16 HHMNSS.S, hours, minutes, seconds, GMT 

17 time, seconds after the hour 

18 number of filters and/or ratios used at this altitude 

19 solar zenith angle, degrees 

20 0.0; for last record, -1. 

Note: There is one type 300. record for each integer altitude 

covered by the Rocoz photometer profile. 


1 400. These records contains ozone content in 12 standard 

pressure layers, derived from the Composite 
Ozone Table given in the type 300. records; where 
no Rocoz data are available, ECC data are used. 

If neither is available, words at that level are 
set to 0.0 

2 1., the first of two type 400. records 

3 air pressure at top of layer, mil 1 iatmospheres 

4 ozone content between above pressure altitude and the 

next pressure altitude, milliatm-cm 
5-20 repeat of words 3 and 4 for the next pressure altitude. 

1 400. Second record, ozone content in layers 

2 2. the second record of this type 

3-8 continuation of (air pressure at top, ozone content) 

9 1000. 

10-20 0.0 

Standard pressure altitudes are at 1000./2(nth power) as follows: 
0.0 0.24 0.49 0.98 1.95 3.9 7.8 15.6 

31. 62. 125. 250. 1000. 
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Word Content 

1 500. Ozone values at 33 standard pressures levels from 

Rocoz Composite Table or, at low altitudes, from 
ECC observations. 

2 level number, 1. through 33. 

3 standard pressure, lowest value for which ozone data are 

available, millibar 

4 altitude corresponding to pressure in preceding word, km 

5 air temperature from datasonde, kelvin 

3 

6 Rocoz Composite ozone, number/m 

7 Rocoz Composite ozone mass mixing ratio 

8 air temperature from ECC sonde, kelvin 

9 ECC sonde ozone, number/m 3 

10 ECC sonde ozone mass mixing ratio 

11 500. 

12 level number 

13-20 repeat of 3-10 for this level, 

or padded with 0.0 to fill out the record 


Standard 

pressure 

levels are. 

in millibars: 

0.05 

0.07 

0.10 

0.15 

0.20 

0.30 

0.40 

0.50 

0.70 

1.0 

1.5 

2.0 

3.0 

4.0 

5.0 

7.0 

10.0 

15.0 

20.0 

30.0 

40.0 

50.0 

70.0 

100.0 

150.0 

200.0 250.0 

300.0 

400.0 

500.0 

700.0 


850.0 1000.0 

If Rocoz, ECC, or datasonde observations are not available at any 
given pressure level, the corresponding words are set to 0.0. 
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